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Reattachment of a Two-Dimensional, 
Incompressible Jet to an Adjacent Flat Plate* 


C. BOURQUE and B. G. NEWMAN?+ 


SUMMARY: As part of a general investigation into Coanda effect, a study has 
been made of the reattachment of a two-dimensional, incompressible, turbulent 
jet to an adjacent, inclined, flat plate. The jet separates from the boundaries 
at the slot lips and reattaches to the plate downstream, a phenomenon which 
is associated with the lowering of the pressure between the jet and the plate 
accompanying the entrainment of fluid there. It is found that the flow becomes 
independent of both the length of the plate and the Reynolds number when 
these parameters are sufficiently large: the flow, scaled with respect to the 
width of the slot, is then uniquely determined by the plate inclination. Two 
approximate theories are developed for the mean pressure within the separation 
bubble, the position of reattachment and the increase in volume flow from the 
slot: the agreement with experiment is fairly satisfactory. These theories are 
a development of Dodds’s analysis for the reattachment of a jet to a plate 
offset from, and parallel to, the axis of the slot and, for the purpose of 
comparison, a limited study is also made of this flow. 


1. Introduction 


Jets are frequently observed to adhere to and flow round nearby solid boun- 
daries. This is seen, for example, when one’s finger is held close to a thin stream 
of water issuing from a tap, or when tea is poured from a badly-designed teapot. 
This general class of phenomena, which may be observed in both liquid and 
gaseous jets, has become known as Coanda effect, after Henri Coanda the inventor, 
who applied it in numerous inventions ’*. In these applications a gas was emitted 
from a thin slot, one wall of which was extended as a series of straight segments 
of increasing length which progressively diverged from the axis of the slot. These 
were arranged so that the flow became reattached to the wall after each discontinuity 
in surface slope and thus the jet was progressively deflected. The deflection of the 


*This report is abstracted from the thesis “Déviation d'un Jet Turbulent Incompressible par un 
Volet Incliné—Effet Coanda™ presented by the first author to the Université Laval for the 
degree of M.Sc. The work was supported by the Defence Research Board of Canada under 
grant number 9550-05. 

+Canadair Visiting Research Professor in Fluid Mechanics, Université Laval, 1958-59. Professor 
of Aerodynamics, McGill University. 
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Incompressible Jet to an Adjacent Flat Plate* 


C. BOURQUE and B. G. NEWMAN* 


SuMMaRY: As part of a general investigation into Coanda effect, a study has 
been made of the reattachment of a two-dimensional, incompressible, turbulent 
jet to an adjacent, inclined, flat plate. The jet separates from the boundaries 
at the slot lips and reattaches to the plate downstream, a phenomenon which 
is associated with the lowering of the pressure between the jet and the plate 
accompanying the entrainment of fluid there. It is found that the flow becomes 
independent of both the length of the plate and the Reynolds number when 
these parameters are sufficiently large: the flow, scaled with respect to the 
width of the slot, is then uniquely determined by the plate inclination. Two 
approximate theories are developed for the mean pressure within the separation 
bubble, the position of reattachment and the increase in volume flow from the 
slot: the agreement with experiment is fairly satisfactory. These theories are 
a development of Dodds’s analysis for the reattachment of a jet to a plate 
offset from, and parallel to, the axis of the slot and, for the purpose of 
comparison, a limited study is also made of this flow. 


1. Introduction 


Jets are frequently observed to adhere to and flow round nearby solid boun- 
daries. This is seen, for example, when one’s finger is held close to a thin stream 
of water issuing from a tap, or when tea is poured from a badly-designed teapot. 
This general class of phenomena, which may be observed in both liquid and 
gaseous jets, has become known as Coanda effect, after Henri Coanda the inventor, 
who applied it in numerous inventions’. In these applications a gas was emitted 
from a thin slot, one wall of which was extended as a series of straight segments 
of increasing length which progressively diverged from the axis of the slot. These 
were arranged so that the flow became reattached to the wall after each discontinuity 
in surface slope and thus the jet was progressively deflected. The deflection of the 


*This report is abstracted from the thesis “Déviation d'un Jet Turbulent Incompressible par un 
Volet Incliné—Effet Coanda” presented by the first author to the Université Laval for the 
degree of M.Sc. The work was supported by the Defence Research Board of Canada under 
grant number 9550-05. 

+Canadair Visiting Research Professor in Fluid Mechanics, Université Laval, 1958-59. Professor 
of Aerodynamics, McGill University. 
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. Approximate mean-flow streamlines for a two-dimensional jet reattaching to an 
inclined flat plate. 


flow at the slot and at each subsequent corner is associated with a reduced wall 
pressure and thus, for given upstream conditions, the mass flow from the slot is 
greater than that of a free jet if the flow is subsonic. Furthermore, the quantity of 
fluid entrained by the jet from the surroundings is apparently increased when the 
jet is deflected. These features were exploited by Coanda in both two-dimensional 
and axisymmetric situations to improve the scavenging of internal combustion 
engines, to augment the thrust of a nozzle with induced flow'*’, to increase the 
maximum lift coefficient of a wing and to achieve possible improvements in the 
performance of other devices. A critical review of these devices has been given by 
Voedisch’. It should be noted that in these applications the surrounding fluid is 
either moving or at rest. 


Although a jet is probably deflected more efficiently by a curved surface'’®’, 
there is nevertheless considerable interest in Coanda’s particular wall geometry, the 
simplest form of which consists of a single flat plate inclined at an angle to the 
axis of a two-dimensional slot. The flow from this slot with the surrounding fluid 
at rest is the main topic of the present paper (Fig. 1). 


After leaving the slot the flow separates from the surface and, if the plate is 
sufficiently long and its inclination is sufficiently small, the jet curves towards the 
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plate and eventually reattaches to if, enclosing a separation bubble. The forces on 
the plate associated with this flow have been measured by von Glahn™ and a 
physical explanation of the behaviour has been given by Squire. Imagine the jet 
flowing freely, i.e. completely separated from the plate. After leaving the slot the 
highly unstable shear layers on both sides of the jet quickly become turbulent and 
the surrounding fluid is rapidly entrained. That fluid which is entrained near the 
plate is slightly accelerated there and thus the static pressure on the wall becomes 
less than that of the surroundings, with the minimum pressure some distance from 
the corner. In consequence the jet curves towards the wall, which further reduces 
the pressure there. Thus eventually, if the wall is long enough, the flow reattaches 
to it, enclosing a separation bubble. It is clear from this explanation that the 
establishment of the flow depends on the maintenance of approximately two- 
dimensional conditions, in particular the prevention of spanwise flow parallel to the 
length of the slot. In practice such flow is prevented by the use of end plates or 
similar devices. It is found experimentally that, at large plate inclinations, the 
flow can be of the type just described or can be fully detached, depending on how 
it is initiated. For example, if reattached flow is established at a low inclination and 
the plate angle is then increased, the flow will remain attached to a relatively high 
angle. On the other hand if the flow is initially detached from the plate and the 
angle of inclination is reduced, the flow will spontaneously reattach to the plate at 
a smaller angle. 


A potential theory for this flow has been given by Metral®’ in which the jet 
is assumed to be fully attached to the wall. The pressure at the sharp corner is 
thus infinitely negative and the theoretical mass flow is too large. The theory is 
developed by examining the flow in the hodograph plane and is similar to analyses 
for flow round other curved surfaces which have been given by Lighthill* and 
Yen”, Of necessity the theory cannot predict the reattachment of the jet to the 
surface or the details of the flow near the sharp corner. In order to do this the 
entrainment of surrounding fluid by the jet must be included in the analysis. 
Entrainment between a shear layer and its surroundings is fundamental to the 
understanding of a variety of viscous flows. For example, an approximate theory 
using the concept of entrainment for a two-dimensional turbulent boundary layer 
has been given by Head“. An entrainment theory for a jet issuing 
parallel to and offset from a flat plate (Fig. 2) has been developed recently by 
Dodds“**. In this latter case the jet reattaches to the plate, enclosing a separation 
bubble, and the theory predicts the mean pressure within the bubble and the 
position of reattachment of the mean flow streamline which divides that flow which 
proceeds downstream along the plate from that which recirculates within the 
bubble. This is the dividing streamline identified by Chapman and his associates 
in their study of separated flows'**’. 


The purpose of the present paper is to apply Dodds’s theory to the Coanda 
deflection of a jet by a flat plate and to compare the predictions with experiment. 


*Note added in proof: A similar analysis has also been made recently by H. R. Chaplin “The 
Effect of Jet Mixing on the Annular Jet”, David Taylor Model Basin Report 1375, 1959. 
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Ficure 2. Approximate mean-flow streamlines for a two-dimensional jet reattaching to an 
offset parallel plate. 


In doing this the original theory has been extended and, for the purpose of compari- 
son, some measurements have also been made of the flow which Dodds studied. 
In addition both flows have been examined from the viewpoint of dimensional 
analysis to isolate the relative importance of the various parameters. The experi- 
mental investigation has been confined to an air jet at speeds low enough for the 
flow to be effectively incompressible. 


NOTATION 
b_ width of the slot 


C, pressure coefficient: — Px) / (Po — Pox) 
Crs (ps Px)! (Po Poo) 


distance of the offset plate from the axis of the slot 


jet momentum per unit span of slot (suffixes 1, 2 and 3 refer to the 
momenta shown in Figs. 3 and 4) 


length of the plate 

static pressure 

stagnation pressure of the fluid supplying the jet 
theoretical static pressure within the separation bubble 
static pressure on the surface of the plate 

static pressure of the surrounding fluid at rest 
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Pp 
Po 
Ps 
Ps 
Poo 
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r theoretical radius of the centre line of the reattaching jet 
Q volume flow from the slot with plate in position 
Q, volume flow from the slot for a free jet 
R_ Reynolds number, =[(p, — p,,) 
s arc length measured along the centre line of the jet 
$s, distance of the slot exit from the hypothetical origin of the jet 


t=tanh {ocy/(s+s,)} on the dividing streamline: suffix 1 refers to the 
particular value of y=y, shown in Fig. 4 


u mean velocity parallel to the axis of the jet or parallel to the plate after 
reattachment 


Um maximum value of wu in the wall jet after reattachment 
x distance measured along the plate from the corner in each case 
Xz Value of x at the reattachment point 
y distance from the jet axis 
z_ distance from the plate 
Zm/2 larger value of z for which u=u,,/2 in the reattached wall jet 
Zm VWalue of z at u=Um 
a angle between the inclined plate and the axis of the slot 


6 angle at which the centre line of the jet theoretically impinges on the 
plate (Fig. 4) 


vy kinematic viscosity of the fluid 
p density of the fluid 
o parameter for a free turbulent jet (equation (7)) 


2. Dimensional Analysis of the Two Flows 


Consider the flow shown in Fig. 1, in which a two-dimensional jet reattaches 
to an adjacent flat plate inclined at an angle < to the axis of the slot. If the ratio 
of the supply pressure p, to the ambient pressure p,, is less than about 1-1, the flow 
may be treated as incompressible and the following parameters are then sufficient 
to define the flow: — 


P.—Px» the supply pressure measured relative to the surroundings 
b, the width of the slot 
I, the length of the flat plate 
a, the angular inclination of the plate to the axis of the slot 
; p, the density of the fluid 


v, the kinematic viscosity of the fluid. 
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The surface pressure p, at distance x from the sharp corner, when expressed 
non-dimensionally as 


Cyn , 
Po Px 


is a function of four non-dimensional combinations of the parameters and hence 
may be written 


(1) 


Similarly the distance x, of the point of reattachment of the mean dividing 
streamline from the corner may be written 


In some cases it is possible to simplify these equations. If the length of the 
plate / is sufficiently large compared with xp, the surface pressure becomes equal to 
p. before the end of the plate and an arbitrary increase in the length of the plate has 
no effect on the flow near the corner. Thus the parameter //b is no longer impor- 
tant and may be omitted from equations (1) and (2). Furthermore, if the Reynolds 
number R =[(p, — p..) 6? /(pv’)]'’? is sufficiently large, the flow may be expected to 
become insensitive to variations of fluid viscosity, so that R may be omitted from 
the equations"*’. It should be noted that the associated limiting value of R may 
itself be a function of 2. 


With these simplifications, the equations become 


A similar set of equations can be stated for the second flow, in which a two- 
dimensional jet reattaches to a parallel flat plate offset from the jet axis (Fig. 2). 
The parameter D/b replaces the plate inclination 2. 


2 
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As before, when //b and R are sufficiently large, the equations simplify to 


A further simplification is possible when D/b is large. In this case variations 
of the separate parameters b and p, — p,, will not appreciably affect the flow as long 
as their product b(p,—p,,), which is approximately the jet momentum per unit 
length of slot, is maintained constant’ '”’. This is equivalent to stating that the 
separate parameters b and p,—p,, may be replaced by the single parameter 
b(p,—p.) and is only permissible if the slot width is small compared with the 
dimensions of the separation bubble; otherwise, the hypothetical change in slot 
position associated with a change of slot width b would be significant. 


With this simplification, 


When //D and [(p, — p,.) bD/(pv)]'? are each sufficiently large, 


Px) D x 

Xr 


3. Theoretical Analysis of the Two Flows 
3.1. GENERAL 

Two theories are given and both are extensions of Dodds’s analysis'**’ for the 
flow shown in Fig. 2. 


The following assumptions are made: — 
(i) The jet flow is incompressible and the mean flow is two-dimensional. 


(ii) The velocity is uniform at the exit from the slot, i.e. the increase in 
velocity associated with the reduced pressure within the separation bubble is 
neglected. Thus the jet momentum per unit length J=2(p,—p,,) b. 


(iii) The pressure within the separation bubble is uniform and the centre line 
of the jet is a circular arc of radius r. The width of the jet is thus small compared 
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with r, which is similar to Taylor’s assumption for a two-dimensional jet blowing 
across a uniform flow'’”’. 


(iv) The jet entrains a similar amount of fluid from the surroundings as a 
free jet and the velocity distribution is similar to that of a free jet. The jet 
momentum is therefore conserved'?”’. 


(v) Entrainment within the bubble ceases, and the flow divides, where the 
extended centre line of the jet intersects the plate. 


(vi) The force parallel to the plate due to the skin friction of the forward and 
backward flow near reattachment is negligible compared with the jet momentum. 
This is a plausible assumption since the flow tends to become independent of vis- 
cosity at high Reynolds numbers. 


The simplest, and probably the most accurate, analysis of a free jet is that due 
to Gdrtler, who assumed that the turbulent eddy viscosity was constant across the 
jet at each downstream station. 


The distribution of velocity parallel to the axis of the jet is given by"® 


sech S+s,’ 


is the axial distance from the slot, 


is the distance of the slot downstream of the hypothetical origin of 
the jet where the flow emerges with momentum J from an infinitely 
thin slot (b—> 0), 


is distance measured perpendicular to the axis of the jet, 


is a constant. Sufficiently far downstream in a free jet it has the value 
7:7 (Ref. 18). In the present theory, however, ¢ is left as a floating 
constant, which permits some variation in the rate of entrainment due 
to the curvature of the jet. 


The mean volume flow between the centre line of the jet and the reattaching 
streamline is, by definition, constant and equal to 4 (Jb/p)*/?. 


Thus the equation of this streamline is 


3 (s+s,) S+5S, 


(8) 


An approximate value of s, may be obtained by equating the hypothetical total 
volume flow in the jet at s=0 to the efflux from the slot, (Jb/p)'’*. 
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This gives s,=ob/3 and the reattachment streamline is given by 


The two flows are now examined separately. 


3.2. REATTACHMENT OF A JET TO A PLATE PARALLEL TO AND OFFSET FROM 
THE AXIS OF THE SLOT 


The mathematical model of the flow is shown in Fig. 3. 


First Theory 


For the first theory the momentum principle is applied locally at the point of 
reattachment, and the pressure difference between the surroundings and the bubble, 
Poo — Ps, is neglected. 

Thus Jcosé=J,-J,= pu*dy — pu*dy, (10) 

where 4 is the angle at which the jet impinges on the plate, J, is the flow momentum 
proceeding downstream along the plate, J, is the flow momentum returned within 


the bubble at the point of reattachment, and y identifies the dividing streamline and 
is given by equations (8) and (9). 


heh 


REATTACHING STREAMLINE 


FiGurE 3. Mathematical model for an offset parallel plate. 
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Substituting equation (7) in equation (10) gives 


3 1 


where ¢ is given by equation (8). 


The following geometrical relationships can be stated from Fig. 3:— 


Substituting relation (12) in equation (9) gives 


and substituting relation (13) in equation (14) gives 
o(4-1)(1 — COS 0) 
38 


where ¢ is related to 6 by equation (11). Hence a unique value of b/r may be 
determined for each value of D/b. 


b 36 


D 


The pressure difference across the jet is, by assumption (iii) of Section 3.1, 


Px 
and thus 
Po — Px r 


The distance of the point of reattachment of the dividing streamline from the 
corner is given approximately by 


1 
y (3 = 1) sin tanh~' 


30 ~ 3f sin 6 ’ 


where ¢ and @ are known in terms of D/b from equations (11) and (15). 


Thus, both C,, and x,/b may be determined as functions of D/b. The in- 
creased volume flow from the slot, which was neglected as one of the original 
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assumptions, may nevertheless be computed approximately. For small values of 
Cys the ratio of the volume flow Q to the volume flow Q, of a free jet with the same 
values of p, — p,, and b is 


4c, 


Second Theory 


For the second theory the momentum principle is applied to the control 
volume WXYZ (Fig. 3) and in this case account is taken of the force arising from 
the pressure difference p,,—ps. It is assumed that the bubble pressure pz, acts to 
the centre line of the jet and then changes discontinuously to the pressure of the 
surroundings, p. 


The momentum equation is therefore 
b 
5). 


ps before the end of the plate and an arbitrary increase in the length of the plate has 
Introducing p,,—Ps=J/r and equation (7), this becomes 


git cos 6, 
from equation (13). 
| 


Equation (19) replaces equation (11) in the first theory: the remaining 
equations are carried over without change. 


Once again Cys, X2/b and Q/Q, may be determined as functions of D/b. 


It should be noted that both the theories may be expected to be inaccurate at 
small values of D/b; this is due, in particular, to the failure of assumptions (iii) and 
(v) of Section 3.1. 


3.3. THE REATTACHMENT OF A JET TO AN INCLINED FLAT PLATE 
The mathematical model for this flow is shown in Fig. 4. 


First Theory 

The momentum principle is once again applied locally at the point of reattach- 
ment and equation (11), giving the angle @ at which the jet impinges on the plate, 
is unaltered. 


The geometrical relationships are different, however. Referring to Fig. 4, the 
arc OA =s=r (z+), where 2 is the angle of inclination of the plate. It is possible 
to proceed by stating the exact geometrical relationship between 4 and 2; however, 
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Ficure 4. Mathematical model for inclined plate. 


since the assumptions of the analysis are such that it is only valid if r is very much 
greater than the jet width, it is reasonable to put @=a. The small error introduced 
fortunately reduces the error associated with assumption (v) and results in a more 
accurate theory for 2 > 20°. 


Substituting this in equations (9) and (16), 


(5-1) 
where ¢ is related to 2 by equation (11), which becomes 


This distance of the point of reattachment from the sharp corner is given approxi- 
mately by 


Xn € tanh-'t 


b 3a sin a (23) 
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The increased volume flow, as in the previous case, is 


Second Theory 


This is similar to the second theory for the previous flow, the momentum 
principle being applied to the control volume WXYZ shown in Fig. 4. Once again 
the bubble pressure is assumed to act to the centre line of the jet and there to 
change discontinuously to p,.. 


In applying the momentum theorem, account must now be taken of the 
momentum J, flowing back along the plate. 


Ps) r(l-—cos2) . (24) 


Since the pressure within the bubble is constant, J, equals the momentum 
inboard of the reattaching or dividing streamline at the plane WX. 


oo 


where y, identifies the position of the dividing streamline as it crosses plane WX. 


With the aid of equations (7) and (25), and introducing p,,— ps=J/r, equation 
(24) reduces to 


avy, 


1 3 3 = 
cos a= where ¢,=tanb ra+S, 


(26) 


Since the arc distance OA is twice the arc distance OB, equation (9) gives the 
following relation between ¢ and f,:— 


(27) 


Equations (26) and (27) replace equation (22) of the first theory and the re- 
maining expressions for C,s, x2/b and Q/Q, (equations (21), (23) and (18)) are 
unchanged. 


Once again it should be noted that both these theories may be expected to 
become inaccurate at small values of 2, for then the centre line of the jet is no 
longer approximately circular, the pressure difference p,.— ps differs significantly 
from J/r and the entrainment within the bubble ceases relatively close to the 
corner. 
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FicurE 5. Reattachment distance for the offset parallel plate. /=12 in. 


4. Details of the Experimental Investigation 


The apparatus which was used to investigate the two flows is described in the 
Appendix. For the inclined plate the investigation was comprehensive and inclu- 
ded measurements of the surface pressure p,, the reattachment distance xz and the 
increase in exit volume flow Q/Q, for various slot widths b, plate lengths /, plate 
angles 2, and supply pressures p,—p.. In addition the range of values of 2 and 
1/b for which either fully-separated or reattached flow is possible, was determined 
at high Reynolds numbers. The velocity distribution downstream of reattachment 
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FiGurE 6. Alternative presentation of the reattachment distance for the offset parallel plate. 
in. 


was also measured. For the offset parallel plate the surface pressure and reattach- 
ment distance were measured for various slot widths, offset distances D and supply 


pressures, using a “long” plate. The effect of plate length was not investigated as 
the main purpose of these measurements was to obtain a comparison with the 
theory valid for long plates and high Reynolds numbers. Undoubtedly for this 
case also there exists a range of D/b and //b for which both fully-separated and 
reattached flows are possible, but this was not investigated. 


All the measurements were made at incompressible speeds with p,/ p< 1-05. 


5. Discussion of the Experimental Results and Comparison with 
Theory 


5.1. REATTACHMENT OF A JET TO A PARALLEL FLAT PLATE OFFSET FROM THE AXIS 
OF THE JeT (Fig. 2) 


Figure 5 gives the distance from the slot to the point of reattachment of the 
dividing streamline for various slot widths b, offset distances D, and supply pres- 
sures (p,—p.). The results have been plotted non-dimensionally in accordance 
with equation (4) and indicate that for these measurements the flow was indepen- 
dent of both //b and the Reynolds numbers [(p, —p,,) b?/(pv*)]'/*, when the latter 
were sufficiently large, in agreement with equation (4a). The limiting value of //b, 
which is a function of D/b, has not been determined. All that can be said is that 
the plate length is unimportant when / is greater than 3-8 xp, but this is unlikely to 
be the lower limit. The limiting value of R varies from about 6,000 at D/b=12 
to 2,500 at D/b=S0. 
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Figure 7. Surface pressure distribution at high Reynolds numbers for the offset parallel 
plate. /=12in. 
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Ficure 8. Alternative presentation of the surface pressure distribution on the offset parallel 
plate. /=12 in. 
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o/d=16-2 

212-01N. 9450 
b = O12 

FIGURE 9. Approximate contours of constant pressure within the separation bubble. 


For large values of D/b, the dimensional analysis indicates that x,/D becomes 
a constant at large //D and [(p, — p,,) bD/(pv*)]''? (equation (6a)). The results are 
replotted in these co-ordinates in Fig. 6 and it appears that the limiting value of 
D/b is about 35 and the limiting value of the appropriate Reynolds number 
[(p, — Pp) bD/(pv*)}'/? is about 2x 10* and appears to be approximately constant 


even for small values of D/b. It should be noted that the fact that these results 
plot in accordance with the dimensional analysis is an indication that the flow was 
effectively two-dimensional, for otherwise an additional parameter, slot length 
divided by slot width, would be of importance. 


The distribution of pressure on the surface of the plate, plotted as C,,= 
(ps — Pp») /(Po — Px») is given in Fig. 7 for various values of D/b when //b and R are 
large. Once again the results are in accordance with the simplified dimensional 
analysis (equation (3a)). Fig. 8 also indicates that when D/b exceeds about 35 the 
pressure distribution plotted in the form {(p,—p)D/(p,—p)6} is a function of 
x/D only and is independent of D/b in accordance with equation (5a). 


The contours of constant static pressure within the separation bubble for 
D/b=16-2 are shown in Fig. 9. The pressure was measured by means of a plate 
fitted with pressure taps and mounted in the centre of the jet in the plane of two- 
dimensional flow (see Appendix). No great accuracy can be claimed for these 
measurements: the pressure taps were of necessity relatively close to the leading 
edge of the plate and furthermore it appears that the secondary flow associated with 
the presence of the plate caused the length of the bubble to decrease by about 15 
per cent. Nevertheless the general shape of the contour is no doubt a good indica- 
tion of the true pressure distribution. The estimated position of the jet and the 
reattaching streamline are indicated in dotted lines in Fig. 9 and are based on the 
known position of reattachment. The mean streamlines within the bubble are 


August 1960 217 


+0086 
\ 0-071 
* +0043 
Cy +0°028 
4 
+0-014 
/ 
/ 
O086 4 
4 4 
4 
-0-0 
=0-05 


Cc BOURQUE AND B. G NEWMAN 


RST THEORY 
SECOND THEDRY 


eibigie 


FicurE 10. Comparison between the theoretical bubble pressure and the minimum pressure on 
the offset parallel plate. /=12 in. 


indicated qualitatively in Fig. 2 and are similar to those associated with a vortex 
situated near the point of reattachment. The back flow within the bubble 
decelerates on approaching the rear wall which joins the plate to the slot, and thus 
the pressure rises into the corner. The pressure on the surface of the plate is 
smallest opposite the centre of the hypothetical vortex. It is seen that this physical 


model of the flow satisfactorily explains the pressure distribution within the bubble 
shown in Fig. 9. It should be noted that the point of reattachment only corresponds 
to the position of maximum surface pressure when the jet impinges on the plate 
at an angle approaching 90° (i.e. D/b is large). For small D/b it more nearly 
corresponds to the position where the pressure coefficient changes sign. 


In relating the bubble pressure and reattachment distance to the theories in 
Section 3.2, it is evident that the basic assumption of a uniform pressure within 
the bubble in incorrect. However the minimum pressure on this surface is roughly 
the mean pressure within the bubble (Fig. 9) and therefore this particular pressure 
is chosen for comparison with the theories in Fig. 10. The measured reattachment 
distances are also compared in Fig. 5. It is apparent that the first theory, with o 
having the plausible value of 12, compares very favourably with the reattachment 
distance, and is approximately correct for the minimum surface pressure, although 
smaller values of 7 give better agreement for the middle range of D/b. The second 
theory is less satisfactory and requires very low values of oc: this is 
associated with the assumption that the low bubble pressure acts on the rear wall 
which joins the plate to the slot, whereas in reality the pressure on this wall is much 
higher than the mean bubble pressure*. Both theories are incorrect when D/b —> 0, 
for then the minimum C,, ~ 0 and is not a large negative value as predicted. 


*Note added in proof: Indeed Chaplin has effectively combined the two theories to predict 
the pressure on the rear wall. 
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Applying the momentum principle perpendicular to the plate suggests that 


uo 


| (x/b)=0; 


however the measured distribution gives non-zero values for this integral. This is 
thought to be associated with the momentum of the entrained fluid perpendicular 
to the plate and to a much lesser extent to the skin friction on the rear wall and 
the outgoing momentum perpendicular to the plate in the wall jet. It is interesting 
to note that Sigalla'*’ also detected a moderately high pressure near the slot for a 
plane wall jet, presumably for the same reason. It should be noted that the neglect 
of this momentum has no effect on the first theory and a small, but probably 
negligible, effect on the second theory. 


5.2. REATTACHMENT OF A JET TO AN INCLINED PLATE 

As mentioned in the Introduction, there is a range of values of 2, depending 
on |/b, for which the jet may or may not be reattached to the plate. The angle 
at which a free jet spontaneously reattaches to the plate, and the maximum angle 
for which reattached flow is possible, are shown in Fig. 11. These measurements 
were made at a sufficiently high Reynolds number for the flow to be insensitive 
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Ficure 11. Zones of flow 
at high Reynolds numbers 
for the inclined plate. 
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to variations in this parameter. The two curves shown divide the figure into three 
regions. To the left the jet is always reattached, forming a separation bubble, and 
to the right it is always separated as a free jet. In the central region, the jet is 
separated or reattached depending on how the flow is started. For example, 
reattached flow may always be obtained in this region by guiding the flow round 
the sharp corner temporarily with an additional surface: conversely it can always 
be separated by temporarily introducing a large disturbance downstream. For all 
the remaining work which is to be described, the flow in the central region is 
reattached. 


The reattachment distances x, for various slot widths b, plate lengths / and 
supply pressures (p, — p,.) are shown in Fig. 12 as a function of the angle of inclina- 
tion of the plate z. The results are presented in terms of the non-dimensional 
parameters given in equation (2). Only the results for the lowest Reynolds number 
deviate significantly from the other results, indicating that R is no longer a 
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significant parameter when it is sufficiently large. Furthermore all the results lie on 
one curve (i.e. are insensitive to //b) when x, is smaller than about 0-7/ for all values 
of z and thus the limiting value of //b at any given a may be detemined as 1-4 times 
the value of x,/b given by this curve. The results are therefore in agreement with 
the prediction of equation (2a). No results were obtained for 2 < 30° because the 
bubble of separation was too small, even with the largest available slot width. 


The distributions of pressure on the surface of the plate for z=35°, 45°, 55°, 
624° are shown in Figs. 13, 14, 15, 16. The first three figures are for values of 
I/b sufficiently large for this parameter to be unimportant and thus show the effect 
of changes of R only. It is noted that when R exceeds about 5,500 the curves 
become independent of R (equation (la)). Once again the results indicate that the 
flow was effectively two-dimensional. With one exception the results shown in 
Fig. 16 were all made at the same high R and thus show the effect of variations of 


Ficure 13. Surface pres- 

sure distribution on the 

inclined plate. /=12 in., 
a= 35°. 
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Figure 14. Surface pressure distribution on the inclined plate. /=12 in., 2=45°. 

| 


12 20 24 28 
a/b 


Ficure 15. Surface pressure distribution on the inclined plate. /=12 in., z=55°. 


l/b only. Reference to Fig. 12 indicates that at 2=624° the flow is certainly 
influenced by the length of the plate when //b is less than 200. 


The universal reattachment distance for large R and //b is compared with the 
two theories for various values of o in Fig. 17. As all the theoretical curves cross 
the experimental curve, it is apparent that none of them strictly predicts the re- 
attachment distance at all angles for one value of 7. Although the appropriate 
values of o will not be identical with that for a free jet, because of the flow curva- 
ture and possibly altered transition in the jet, nevertheless they are not expected to 
differ appreciably from the value 7.7. On this basis the second theory is preferable 
to the first. 
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FiGuRE 16. Surface pressure distribution on the inclined plate. /=12 in., 
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FiGurRE 17. The reattachment distance on the inclined plate: comparison with theory. 
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FIGURE 18. Pressure within the separation bubble for the inclined plate. /=12 in. 


The variation of surface pressure in the neighbourhood of the slot is smaller 
than in the previous case and thus the assumption of a constant bubble pressure is 
more acceptable. The point of reattachment corresponds to the position of maxi- 
mum surface pressure for large z and to the position where C,=0 for small values 
of 2, which behaviour is similar to the first flow, and the general distribution of 
pressure may be associated with the curvature of the streamlines as before. 
Unfortunately it was impossible to measure the distribution of pressure within the 
bubble because, when the bubble was sufficiently large, the pressure differences 
were small and highly oscillatory. For comparison with the two theories, both the 
minimum pressure on the surface and the pressure at the slot determined by extra- 
polation, have been plotted in Fig. 18. Once again the second theory reauires 
values of « closer to those of a free jet. 


Measurements of the increase in volume flow due to the presence of the plate 
are given as a function of z in Fig. 19. The increase is very small and could not 
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be measured with great accuracy. It is interesting to note that the maximum 
Two theoretical curves from the second theory, for 


increase occurs at 2 =~ 25°. 
the values of « which gave best agreement with the pressure at the corner, are also 
shown in the figure and the agreement is fair for values of z between 20° and 50°. 
The theory is in error at low values of 2 because the basic assumptions fail. It is 
also possible to compare the results with the potential flow analysis due to Metral? 
and it is noted that this theory gives increments which are far in excess of those 
Viewed in this light, the present entrainment theory is a 


achieved in practice. 
decided improvement for 2 greater than about 10°. 
In this case the application of the principle of momentum perpendicular to 


the plate suggests that, for large //b, 


d (x/b)=2 sin z. 


0 


However, the integrated surface pressure distributions are, in general, less than 
2 sin 2, although the agreement is fairly good for the largest value of z. Once again 
therefore the incoming momentum of the entrained fluid perpendicular to the jet 


is apparently significant. 
Referring to Figs. 12, 13, 14 and 15, the effect of a reduction of R below those 
values for which it is an insignificant parameter, is to reduce the reattachment 
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distance x,/b and to increase the numerical value of the minimum C,, within the 
bubble, both of which are associated with an increased curvature of the jet. The 
reason for this variation is not well understood. It is likely to be associated with 
the skin friction in the region of reattachment, but a variation in the right direction 
would be obtained only if the total skin friction associated with the back flow was 
greater than that associated with the flow proceeding downstream in the vicinity of 
the reattachment point. 


In Figs. 12 and 16 it is seen that, for a given value of 2, the effect of a reduction 
of //b below the limiting value is to reduce x,/b and to increase (—C,,.) in the 
bubble: it is thus similar to the effect of a reduction of R. In this case the 
explanation is as follows. When the point of reattachment is near the end of the 
plate the final direction of the flow is inclined at an angle greater than 2 to the axis 
of the slot, because the flow outboard of the reattaching streamline is no longer 
constrained to depart in the direction of the plate. This statement is supported by 
the fact that 


[ Cy. d (x/b) 


0 


in Fig. 16 increases as //b decreases. Thus the final momentum in the direction 
of the plate is reduced. According to the theoretical ideas of Section 3.3, this is 
associated with a reduction of both the bubble pressure and the reattachment 
distance. 


It was observed that the distribution of surface pressure and of the position 
of reattachment were fairly steady in time for 2 << 35°, but fluctuated with increas- 
ing amplitude as 2 was increased, until the reattachment point reached the limiting 
value 0-7/. From then on the fluctuations decreased until the flow finally failed to 
reattach to the plate. The magnitude of the maximum fluctuations in xz was +20 
per cent. The oscillations are due to the high turbulence in the jet and to a lesser 
extent to variations in the position of transition, both of which lead to fluctuations 
in the effective value of the plate angle z. The fluctuations are therefore roughly 
proportional to dx,/d2 and it is interesting to note that, on this basis, the descrip- 
tion just given of the increase and decrease of the fluctuations with 2 is in agree- 
ment with the curves for each value of //b shown in Fig. 12. 


Referring to the right hand curve in Fig. 11, it is noted that the maximum 
possible plate angle for reattached flow does not vary a great deal with //b. For 
example, for //b=25, amax=57°, while for 1/b=200, 2msx=65°. The maximum 
angle which has so far been obtained experimentally is 75°, corresponding to an 
[/b of approximately 2,000. The theories of course predict 90° for an infinitely 
long plate. 


Pitot traverses of the reattached flow have been made at angles 2=30° and 
45° for | xg. The traverses were made far enough downstream for the static 
pressure to be sensibly atmospheric, as indicated by Figs. 13 and 14. The 
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FiGurE 20. Velocity profiles downstream of reattachment. 


computed velocity profiles are presented non-dimensionally in Fig. 20 in terms of the 
maximum velocity u,, and the distance Z,,., which is the larger value of z at which 
u=4u,. Plotted in this way the profiles are similar and in good agreement with 
the mean profile representing the measurements of a wall jet by Forthmann'’*’ and 
by Sigalla®® at slightly higher Reynolds numbers and with Glauert’s'*’) theory for 
Zm/Zm/2=0°14, where Z,, is the value of z at u=u,. Thus it appears that a conven- 
tional wall jet is obtained far downstream of reattachment. 


6. Conclusions 


The reattachment of an incompressible, two-dimensional jet to a neighbouring 
surface is associated with the entrainment of surrounding fluid by the jet. For a jet 
reattaching to an adjacent inclined flat plate, it is found that the flow becomes 
independent of the Reynolds number at the slot, [(p, — p,.)b?/(pv’)]'/*, when this 
exceeds about 5,500 and independent of the non-dimensional length of the plate, 
I/b, when this exceeds a value which depends on the plate angle « and which may 
be determined from Fig. 12. For a jet reattaching to an offset parallel flat plate 
the flow becomes independent of D/b, !/b and Reynolds number when these are 
sufficiently large, the limiting value of D/b being about 35. 
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Approximate theories have been developed for the mean pressure within the 
separation bubble and the position of reattachment. These are based on the known 
properties of a free jet and the momentum principle applied to a control volume. 
For the point of reattachment, the agreement with experiment is very satisfactory 
for the offset plate, using the first theory, for D/b > 3 and is fairly satisfactory for 
the inclined plate, using the second theory, for 2 > 25°. 


The increase in volume flow from the slot associated with the low pressure 
on the inclined plate has a maximum value of 5 per cent for a long plate and occurs 
at a= 25°. This is very much less than that predicted by potential flow theory. 


The flows exhibit a hysteresis phenomenon. For the inclined plate there is a 
range of 2 at each //b for which the jet may be either reattached or fully separated, 
depending on how the flow is started (Fig. 11). For the offset plate the same is 
presumably true with D/b replacing 2, although this was not investigated. The 
phenomenon may be of importance when deciding on the plate length in any 
practical application. 


The present theories could usefully be improved, particularly for small values 
of z or D/b, and could also be extended to the case of a compressible subsonic jet. 
It might also be interesting to examine the importance of entrainment for the 
deflection of a sonic or a supersonic jet round a sharp corner, although the existing 
potential flow theories’: *’ are probably more satisfactory and there is, of course, 
no increase in mass flow in this case. 
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FiGcurE 21. Sketch of the apparatus. 
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Ficure 22. Velo- 
city profile at the 
exit from the slot 
with a free jet. 


Appendix 


APPARATUS AND EXPERIMENTAL PROCEDURE 


The apparatus, shown schematically in Fig. 21, was designed to produce a quasi 
two-dimensional jet under those conditions which have been examined theoretically. 
Air was supplied from a centrifugal compressor driven by a constant speed, three-phase 
motor. The flow to the apparatus was controlled by a bleed valve situated upstream of a 
large reservoir: the reservoir reduced the fluctuations in total pressure to an acceptable 
level of less than one per cent. 


In order to obtain a uniform velocity at the slot, the flow was decelerated by means 
of a 40° expansion, in which a perforated plate was placed to inhibit separation of the 
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flow from the walls. This was followed by a short uniform settling space which contained 
a deep-cell honeycomb. The flow was then accelerated into the slot, the width of which 
was adjustable: the minimum value of the contraction ratio, corresponding to the 
maximum slot width, was 11:1. This arrangement proved to be very satisfactory; the 
pressures given by three pressure tappings across the settling space were identical and 
equal to the total pressure in the centre of the jet, which was itself uniform along the 
length of the slot. With the flow issuing as a free jet, Pitot traverses at the slot indicated 
that the flow emerged with laminar boundary layers about 0-010 in. thick. Traverses for 
two slot widths are shown in Fig. 22. 


The effective length of the slot was 21 in. and its width could be varied up to a 
maximum value of 0-19 in. The lips were surface ground and the slot width in the centre 
was less than that at the ends by 0-005 in. with wind off and, at the highest stagnation 
pressure (20 in. water), was greater by 0-003 in., due mainly to deformation of the top lip. 
Most of the readings were made at an intermediate stagnation pressure and the variation 
of slot width was roughly 0-002 in. 


Three plates were used, respectively 3, 6 and 12 in. long. They were cut from 
“Perspex” plastic sheet, } in. thick. Pressure tappings of 0-015 in. diameter were drilled 
in the longest plate, the holes being closest together near the slot. The plywood sides of 
the setting space were extended to form large end walls. These were insufficient, however, 
to completely prevent inflow from the sides, as detected by tufts, and it was necessary to 
attach additional end walls directly to the plates in order to produce satisfactory two- 
dimensional conditions. 


Figure 23. Mounting for a 
parallel offset plate. 


The plates could be mounted either parallel to the axis of the slot and some distance 
from it (Figs. 2 and 23) or inclined to the axis to form a sharp corner with the top lip of 
the slot (Figs. | and 21). The outside of the top lip and the leading edge of the plates were 
chamfered at 45° so that they could abut to form the sharp corner in the latter case. 
The leading edge of the plate was held in position against the suction force on the plate 
by means of the four nylon strings shown in Fig. 21. The sharp corner was sealed at 
the back of the plate with adhesive tape. 


In each case a desired slot width was set with feeler gauges. The pressure in the 
settling space was measured on a vertical alcohol manometer to an accuracy of one 
per cent. The surface pressure on the plate was measured with a vertical multi-tube 
manometer, filled with alcohol and fitted with a freezing device. The oscillation of those 
tubes connected to the tappings near reattachment were about 20 per cent of the mean, 
and care was taken to freeze the manometer to give approximately the mean pressure 
distribution. 
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Ficure 24. Paper flags for determining the position of reattachment. 


The position of reattachment of the dividing streamline was determined by using 
small flexible paper flags (two-dimensional tufts) fixed to the surface with hinge lines 
perpendicular to the direction of mean flow (Fig. 24). These indicated the instantaneous 
division of forward and backward flow and the mean position of reattachment could 
readily be estimated. 


The mass flow through the slot was measured in terms of the pressure difference 
between a central Pitot tube and a wall static tapping in the circular pipe which connected 
the apparatus to the large reservoir. 


The velocity profile downstream of reattachment was measured with a round 
Pitot tube with internally sharpened lips and external diameter 0-028 in. This was 
mounted on a dial-gauge traversing gear. The traverses were made where the static 
pressure was sensibly equal to the pressure of the surroundings. 


The pressure distribution within the separation bubble for the offset plate was 
determined with an additional } in. thick plate fitted with three pressure tappings and 
rounded leading edge, mounted in the centre of the plate and parallel to the end walls. 
The taps were situated at 4 in. from the leading edge. 
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Analysis of Stress Concentrations at Reinforced 
Holes in Infinite Sheets 


W. H. WITTRICK, M.A., Ph.D., F.A.A., A.F.R.Ae.S., A.F.LA.S. 


(Department of Aeronautical Engineering, University of Sydney) 


SUMMARY: Equations are derived which enable the stresses around a compactly 
reinforced hole in an infinite sheet to be computed. The only restrictions 
placed on the shape of the hole are that it has at least one axis of symmetry 
and that the region outside it can be conformally mapped on to the region 
outside the unit circle by a transformation function in the form of a poly- 
nomial of any order. This admits a very wide range of holes of practical 
importance, including circles, ellipses, and squares, triangles and rectangles 
with rounded corners. Two loading cases are treated; the first corresponds to 
uniform tensions at infinity in directions parallel and perpendicular to the axis 
of symmetry of the hole, and the second to uniform shear at infinity in these 
directions. Superposition of these basic loading cases enables the stresses to be 
determined for any uniform state of stress at infinity. The equations are in 
an ideal form for use with a digital computer. 


Introduction 


Considerable interest attaches to the problem of determining the stress 
concentration arising at a hole in an infinite sheet, subjected to stresses at infinity. 
In practice a reinforcing member is usually provided around the hole boundary in 
order to reduce the stress concentration. 


Mansfield’ has shown how, for any given state of stress in the sheet, the hole 
may be designed in such a way that its presence leaves the state of stress unchanged. 
The hole is then said to be a neutral one, and its shape and the distribution of 
reinforcement around its boundary are dictated by the state of stress for which it is 
designed. With any other loading the hole is no longer neutral. Mansfield’s theory 
leads in some cases to practically possible holes, the best known example being for 
a uniform state of stress with both principal stresses of the same sign. The hole is 
then elliptical, with a ratio of axes equal to the square root of the ratio of the 
principal stresses and the major axis in the direction of the greater principal stress. 
Such a hole could be used, for example, for a window in a pressure cabin, in which 
the predominant loading is due to the pressurisation, with a hoop tensile stress 
roughly equal to twice the longitudinal stress. 
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In practice, however, neutral holes suffer from three disadvantages. First, the 
reinforcement required is invariably non-uniform and therefore undesirable from a 
manufacturing point of view. Second, the weight of reinforcement required is often 
rather excessive and it is known that a smaller amount of uniform reinforcement 
can lead to holes having stress concentration factors only slightly in excess of unity. 
Finally, the shape of the neutral hole may be unfavourable for reasons quite apart 
from strength considerations. There is little doubt that these arguments, and 
particularly the weight penalty which must be paid for a neutral hole, have led 
designers to accept compromise solutions, with different shapes and a uniform 
reinforcement having as small a weight as possible consistent with an acceptable 
stress concentration factor. 


At the present time only elliptical and circular reinforced holes have been 
analysed to any great extent® * and the information available to the designer for 
other shapes of hole is rather meagre. The analysis of this paper is intended as a 
means of rectifying this deficiency. It has been made as general as possible, to cover 
holes of different shapes, with any reinforcement, and with any uniform state of 
stress at infinity. 


The only assumptions made about the shape of the hole are, first, that it has at 
least one axis of symmetry and, second, that the region outside the hole may be 
mapped on to the region outside the unit circle by means of a polynomial trans- 
formation function of arbitrary order. The assumption of symmetry is not a serious 
restriction, since the great majority of holes of interest in practice have the required 
amount of symmetry. The second assumption is more restrictive but, nevertheless, 
a large number of practical shapes can be approximated to by this type of 
transformation function, particularly squares and triangles with rounded corners“. 


The reinforcement is assumed to be compact, which means that its flexural 
rigidity is sufficiently small for the shear force in the reinforcement to be neglected 
in comparison with the tension. This forms the basis for most previous work on 
reinforced holes. 


The stress system at infinity is assumed to be uniform but otherwise arbitrary. 
However, it was found to be convenient to analyse separately the two cases of 
tensions at infinity, parallel and perpendicular to the axis of symmetry of the hole, 
and shear at infinity in these directions. The principle of superposition then enables 
any state of stress at infinity to be analysed. 


The equations given are in an ideal form for solving on a digital computer. 
In fact it is possible to write a general programme which will enable the stresses to 
be computed around an arbitrary hole, conforming to the given restrictions, merely 
by specifying the coefficients of the polynomial transformation function, the cross- 
sectional area of the reinforcement, and the state of stress at infinity. 


NOTATION 
x,y Cartesian co-ordinates 


$,n co-ordinates tangential and normal to the hole boundary 
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z=x+iy 

t thickness of the sheet 
A cross-sectional area of the reinforcement 

r local radius of curvature of the hole boundary 


X,Y co-ordinates of points on the hole boundary, made non-dimensional 
by dividing by R 


G_ shear modulus 
~©Poisson’s ratio 
€ angle between the x-axis and the tangent to the hole 
{ complex co-ordinate in the circle plane 
6 polar angle in the circle plane 
ao =e", complex co-ordinate on the unit circle 
transformation function 
R_ constant determining the size of the hole, see equation (6) 


c; real constants in the transformation function (j= —1, 0, 1,.... 
See equation (6) 


F Airy stress function 


?, (2), Xo (Z) complex stress functions, see equation (1) 
(z)= x," (z) 
(2) 
(Q=W, (2) 


a,,6, real constants in the expansions of ‘< ({) and WV (0), for the case of 
tensions at infinity, see equations (25) and (26) 


Px Gn Teal constants in the expansions of ?({) and WV (Q), for the case of 
shear at infinity, see equations (39) and (40) 


f,.f. applied tensile stresses at infinity in the x, vy directions 
S applied shear stress at infinity in the x, y directions 
T tension in the reinforcement 
P= X?+Y? 


+(1 


A prime (’) denotes differentiation 
A dot denotes differentiation with respect to 6 


A bar denotes the complex conjugate value. 
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2. Equations of Plane Stress in Complex Form 


In this section the basic equations of the two-dimensional theory of elasticity 
in complex form“? will be summarised. 


In terms of the complex variables z and z, defined by the equations 
z=x+iy, z=x-iy, 


the biharmonic equation for the Airy stress function F becomes 


Its general (real) solution is 
F=Z®, (2)+zP, (2+ W 


where ‘P, (z) and x, (z) are two arbitrary analytic functions, and a bar denotes the 
complex conjugate. 


The stresses x, fp and XP are given by the equations 


+2 
=4 =4[®,’ (z)+,’(z)],. ‘ (2) 


where WV, (z)=x, (2) 
and a prime (’) denotes differentiation. 


The displacements u and wv of a point in the x and y directions are given by 
the equation 


G (u+ iv)=x®, (z)-z®,’ (z)-V, (z), (4) 


where G is the shear modulus and « is a constant depending on the value of 
Poisson’s ratio ». For the case of plane stress 


3-u 
, 
The only other basic relation required expresses the fact that the sum of the 


direct stresses at a point in any two perpendicular directions is constant, independent 
of the orientation; thus 


-_ 


33+ An =IR+ Fp, (5) 


where s and nm are any two perpendicular directions. 
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3. Transformation of the Equations 


Let the region outside the hole in the z-plane be mapped on to the region 
outside the unit circle in the (-plane by means of a transformation function 


‘ ‘ ‘ (6) 


in which R is a constant, having the dimensions of a length, whose magnitude 
determines the size of the hole, and the dimensionless coefficients c, determine the 
shape of the hole. As mentioned previously, it will be assumed that the hole has 
at least one axis of symmetry, and the x-axis will be chosen to coincide with it. The 
effect of this is to ensure that R and all the c, coefficients are real, thereby simplifying 
the subsequent analysis considerably. In practice and without loss of generality, 
the values c_,= 1 and c,=0 may be postulated. 


Let (2)=, 

and (2=¥, (2 QJ=V 

It can be shown? that the functions @ (¢) and V ({) are expressible in the form 
W(Q=R (BO +b, + (8) 


provided that the stress at infinity is finite and that the resultants in both the x and y 
directions of any tractions which may be applied to the hole boundary are zero. The 
constants z and 8 may be expressed in terms of the stresses at infinity by means of 
the equations 


+h) +iC. ‘ (9) 


(10) 


where f,, f, and S are the values of £X, VY and XV respectively at infinity and C is an 
arbitrary real constant. It can be shown from equation (4) that the values of C, a, 
and b, correspond to rigid body movements which do not affect the stresses and can 
therefore be disregarded. 

It can also be shown? that if the resultant moment, about the origin, of the 
tractions applied to the hole boundary is zero, the real part of the function 


(z)+ z¥, (z) — x, (z) 


\ 


is single-valued. 
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Now the functions ©,’ (z) and Y, (z) are both single-valued, as can be seen from 
equations (7) and (8). This follows because the stress components are single-valued. 
Hence the real part of x, (z) must be single-valued. 


Now Xo (z)= (z) dz= YY’ 


The term in (-' in the integrand of this expression will, when integrated, give 
rise to a term in log ¢ which increases by 27i on making any anti-clockwise circuit 
surrounding the origin. It follows therefore that 


F(b,c_, —Be,)=90, 


where ./ denotes the imaginary part. 


But all the c; coefficients are real and, from equation (10), the imaginary part 
of 8 is 4iS. Hence it follows that 


4. Boundary Conditions 


Consider the equilibrium of an element ds of the boundary reinforcement, with 
a triangular element of the attached sheet, as shown in Fig. 1. It will be assumed 
that the shear force in the reinforcement may be neglected in comparison with the 
tension T. 


Y 4 T+dT 


Ficure 1. Element of reinforcement and attached sheet. 
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For equilibrium of the element in the y direction 


d(T sin £)=t (¥) dx —Xy dy) 


dy) 


where ¢ is the sheet thickness and € is the slope of the boundary. 


Similarly, for equilibrium in the x direction, 
d(T cos £)= —td (@F /éy). 


Multiplication of the first equation by —i and addition to the second gives 
d (Te~*) itd i 2itd 


Integration of this equation and substitution of equation (1) now gives the 
equation 


For future reference it may be noted that the constant of integration has been 
incorporated in the function ®, (z). 


Since T is essentially real, the first boundary condition expresses the fact that 
the right-hand side of equation (12) is real. It will be seen later that this enables 
the function WV, (z) to be determined in terms of ®, (z). 


The secondary boundary condition arises from the equality of strain of the rein- 
forcement and of the sheet to which it is attached, from which the following 
equation is obtained : — 


(13) 


where A is the cross-sectional area of the reinforcement and the s and n directions 
are shown in Fig. 1. Consider now the equilibrium of an element of reinforcement 
alone, as shown in Fig. 2. Resolving in the n and s directions, 


(14) 


(15) 


where r is the local radius of curvature of the boundary. 
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Ficure 2. Element of reinforcement. 


Combining equations (13) and (14), the second boundary condition may now 
be written in the more convenient form 


T[44 |= sst+nmn . ‘ . (16) 
where T is given by equation (12) and (ss + mn) by equations (2) end (5). 


5. Geometry of the Hole 


Some expressions concerning the geometry of the hole will now be derived in 
terms of the transformation function. 


Let the x and y co-ordinates of a point on the boundary of the hole be denoted 
by XR and YR, where X and Y are non-dimensional. 


On the unit circle 
so that, on separating real and imaginary parts of equation (6), 


N N 


X= c,cos j6, Y=— = c,sin j6 
| 1 
Hence X=-— & je,sinjé, Y=- & je,cos jé (19) 
| 1 
N 
and X=-— & cos j6, Y= & fe;sinj. . (20) 
~1 -1 


where dots denote differentiation with respect to 6. 
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The slope é of the boundary of the hole is given by 


(Y do 
= 


do=icd6 and =dé. 


cl (oc) 
(oc)! 


dod (X+iY) 


Also do dé ic 


Hence = RP"?, 


where P=X?+Y?. 


The radius of curvature r of the boundary is readily shown to be given by 
the equation 


6. The Case of Tensions at Infinity 


Consider first the case where the stress system at infinity consists of tensions 
f, and f, in the x and y directions, with no shear stress (S =0). 


From equations (7) to (10), therefore 


It may be noted that although a, does not affect the stresses it is necessary to 
retain it in order to account for the constant of integration involved in the 
derivation of equation (12). 

Because of the symmetry of the hole and of the loading about the x-axis, it is 
clear that if the displacements at a point z, are u, and v,, then at the image point 2, 
the displacements are u, and —v,. Using equation (4) it is then a simple matter to 


show that 
(z)=, (2), (2) =, (2). 


This means that all the a, and b, coefficients in equations (25) and (26) are real. 
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Since §=0 and b, is real, equation (11) is automatically satisfied. 


Now ®,(2)= £4, 


Hence, on using equations (21) and (22), equation (12) for the tension in the 
reinforcement becomes 


TRP'? 


=o (c) [® (7) + + 0D. (0) (0). Q7) 


The right-hand side of this equation must be real. Hence, using the fact 
that  =1/c, the following boundary equation is obtained 


of [ +0 [ (c)+¥ (*)| + 
(8) 
This equation is multiplied throughout by do/[2zi(o—Q], where |{| >1, and 


integrated round the unit circle. Then, on using equations (6), (25) and (26) and the 
known properties of Cauchy integrals, it is found that 


VO=3 (2) ]+ eG )oo- 
1 N-1 N } 


Thus, the function V(¢) has now been determined completely in terms of ® (9), 
epart from the coefficient b, which remains unknown. 


Substitution of equations (6), (25) and (29) into equation (27) now gives the 
following expression for the tension T : — 


N N-1 N 
- a, > jc, cos(n—j) 6| + & [ a, (j-n) c,cos (j—n) 6 | . (30) 
n=0 -j=-1 n=0 j=n+l 
The sum of the direct stresses in the sheet at the boundary of the hole is given by 


H+ fin =4[®,’ (z)+ ®,’ 


(c) 
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Hence, on noting that 
YY (0) Y (6) 
and using equation (22), 


Substitution of equations (6) and (25) now gives the equation 


N N 


- jc; cos (j+ 1) 6+ | na, jc, cos (n—j) 6 | . 
1 1 


Now the strain boundary condition on the hole, equation (16), can be written in 
the form 


1/2 
TP (B+ AM, age 


~2 


Substitution of equations (30) and (31) into (32) then gives 


~ N N-1 
= a, (1 +nQ) = ic, cos (n—j) 6 | + a, (j—n) c, cos (j—n)6 | + 
=0 i 1 n=0 +1 


N 
+ = cos (j+ 1) 4+ (2c cos 26—c,)=0. 


-1 


(34) 


Equation (34) enables the unknowns b, and the a, to be determined. It should 
be noted however that, in general, P, r and therefore Q are functions of 9, which 
makes it impossible to compare coefficients of the various cosine terms directly in 
equation (34). Theoretically Q could be expanded as a cosine series itself and the 
resulting products of cosines in equation (34) could be expressed in terms of the 
sums of cosines. Then, after comparing coefficients, an infinite number of infinite 
simultaneous equations in the unknowns would result. However, even in the simple 
case of an elliptical hole, this process is far too lengthy and complicated to be a 
practical possibility. Furthermore, only an approximate solution could be obtained 
in this way, by retaining only a finite number of a, coefficients and an appropriate 
number of simultaneous equations. 


A far simpler and more direct approach to the problem is to retain only a finite 
number, say m, of the a, coefficients in the first place and to satisfy equation (34) at 
(m+ 1) points in the range 0 < 6< =, say at 9=0, =/m, 2=/m,.. ., =. In this way 
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sufficient simultaneous equations can be derived to determine the values of b, and 
the a,. This procedure is ideal for use with a digital computer and was successfully 
employed for obtaining the stresses around reinforced elliptical holes:*. If it is 
expected that, over some region of the boundary, rapid variations of stress may 
occur, as in a square hole with only slightly rounded corners, the procedure can be 
modified by satisfying equation (34) at smaller intervals of 4 in this region. A good 
guide to whether a sufficient number of the a, coefficients has been retained, for the 
required accuracy, is the convergence of a, to zero as nm increases. Also, if there is 
any doubt about whether the chosen intervals of 4 are small enough, a check can be 
made after the determination of the unknowns to see how closely equation (34) is 
satisfied at the middle points of the intervals. 


There are certain special cases in which Q is constant. For example, if the 
hole is unreinforced, then A=0 and Q=0. If the hole is circular, with a uniform 
reinforcing member, then P, r, A, and therefore Q, are constant. Also, it can be 
shown”? that any elliptical hole designed to be neutral for some particular ratio 
A of f,/f, is of such a shape, and has a reinforcing member with such an area 
distribution, that the variations of P, r, and A with 6 cancel each other out and leave 
Q constant with the value 2/(A+A~'). If, for whatever reason, Q happens to be 
constant, the unknowns can be determined directly from equation (34) by comparing 
coefficients of the various cosine terms. In this case it is easily seen that a,=0 if 
n> N, so that the series for & (0) terminates at a finite number of terms. 


Once the unknown coefficients a, have been determined, the sum of the direct 
stresses, (S$ + 7m), in the sheet at the boundary of the hole may be computed from 
equation (31). The tension 7 in the reinforcement follows from equation (32). 


Thus OP’ 2 (3 + Ani. 


Equation (14) then gives the stress Am and, knowing both this and (53+), the 
stress §§ can be computed. 

Finally, an expression for the shear stress $i on the boundary can be obtained 
from equation (15). 


Now ds = \dz| = (c) do| = RP''*d6. 
l dT 
Hence n= 
_gp-ij2 ‘ mall 
4Rt 
4d 2 TP*!? dP 


P d6\4Rt P* 4Rt dé° 


Using equation (32) this becomes 


dP 
dé 
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Also, from equation (30), 


1/2 N 
G+ DesinG+ 8-5 c_, sin 20+ 


dé \ 4Rt 
N N-1 N 
| a, 6] - | >> (j—n)’c, sin (j—n) | . 
n=0 j= -1 n=0 j=n+l 
(37) 


From equation (23), 


dP 
=2(XX+ YY). 


Equations (36), (37) and (38) enable Ss” to be computed. 


7. The Case of Shear at Infinity 


Consider now the case where the stress system at infinity consists of a pure 
shear stress S, with no tensions (f, =f, =0). 


The symmetry is now such that if the displacements at a point z, are u, and v,, 
then at the image point z, the displacements are —u, and v,. It then follows from 
equation (4) that 

VY, (2). 

Hence the coefficients in equations (7) and (8) are purely imaginary and, on 

using equations (9) and (10), the expansions of <P (Q) and (Q) may be written in 


the form 


where the p, and gq, coefficients are real. 


Comparing equations (8) and (40) it is seen that 


b, =iq,S. 


Hence, equation (11) gives 
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The analysis now follows a course exactly similar to that of the preceding 
section. The fact that the right-hand side of equation (27) is real gives V (0) in 
terms of ®({). Thus 


N-1 N 
n=0 j= 1 


n+ 


The tension 7 in the reinforcement is given by 


4Ri Ssin26+S = Pa je, sin (j n) 6 


N-1 N 
(43) 
n=0 j=n+l1 


The sum of the two direct stresses in the sheet at the boundary of the hole is 
given by the equation 


N 


Substitution of equations (43) and (44) into equation (32) then gives the following 
equation for the determination of the unknowns p,. 


N N-1 N 
| (1 = ic, sin 6 | - [ > (j—n) c,sin (jn) 6 | = : sin 20. 
j= n=0 j=n+l 
(45) 


=0 


Knowing the p, coefficients, the sum of the direct stresses (S§+77) can be 
computed from (44). Equation (35) then gives the tension 7, and 7n follows from 
equation (14). Knowing and nn, the stress S$ can be computed. 


The shear stress 57 at the boundary is given by equation (36) where, in this case, 


d TP"? 
do \4Rt 


) -Scos26+5 iG-me cos (j —n) = 


N-1 N 
-S = > cos (j—n) 6 | . (46) 


n+l 
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As before, since Q is in general a function of @, only an approximate solution is 


possible, by retaining a finite number, say (m — 1), of the p, coefficients and satisfying 
equation (45) at (m-—1) points in the range 0 << 6 <=, say at 2e/m,.. ., 
(m-—1)/m. The equation is automatically satisfied at and =. 


If Q happens to be constant, for any of the reasons previously described, the 


unknowns p, can be obtained from equation (45) by comparing coefficients of the 
various cosine terms. As in the case of tensions at infinity, it is readily seen that 
p.=0 if n> N and the series for P (() terminates at a finite number of terms. 
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Limitations of the Sonic Area Rule 


L. M. SHEPPARD 


(Weapons Research Establishment, South Australia) 


SUMMARY: Recent free flight tests have provided information on the 
applicability of the sonic area rule to non-lifting configurations. An analysis 
of results for the sonic wave drag of swept-back, tapered wings has suggested 
that the sonic area rule is applicable to thin wing and slender body combina- 
tions provided that the product of the wing span/length ratio and the cube root 
of the thickness/chord ratio is less than unity. The wing span/length ratio 
was found to be a much more useful slenderness parameter than the aspect ratio. 


When the wing has a round-nosed section, theoretical considerations predict 
the existence of a so-called leading edge drag force varying with Mach number. 
One would expect such a wing to have the same sonic drag-rise as the equivalent 
body with an identical cross-sectional area distribution. Some test results do 
not confirm this drag-rise equivalence. It appears that the transonic drag-rise of 
the leading edge force is much greater than that predicted by theory. 


Tests on one basic wing-body combination, with different body cross- 
sectional shapes of the same area, suggest that the cross-sectional shape, of a 
smooth slender body, has no influence on the transonic wave drag of a 

configuration. 


1. Introduction 


The sonic area rule’: *’ predicts that the sonic wave drag of a non-lifting wing- 
body combination depends primarily on its cross-sectional area distribution. 
Whitcomb'"’ showed that the rule does enable configurations with low sonic drag-rise 
to be designed. Lord'*’ examined the theoretical design of low drag combinations 
as well as the estimation of sonic wave drag. 


It is very important to know what limitations should be placed on the 
applicability of the sonic area rule, otherwise ineffective designs and incorrect drag 
estimates could occur. Spreiter'*’ and Page investigated the range of applicability 
of the rule for rectangular and delta wings respectively. The range of applicability 
of the sonic area rule for swept-back, tapered wings is investigated here by analysing 
some experimental results obtained by the Aerodynamics Division, Weapons 
Research Establishment’®:*’. Other results obtained at the same time provide 
information on the influence of round-nosed wing sections and the body cross- 


sectional shape. 


The wind tunnel and free flight experimental results utilised in this paper are 
discussed in Section 2. An examination of the range of applicability of the sonic 
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area rule is made in Section 3. Section 4 presents a discussion of some data on wings 
with round-nosed sections. Finally, in Section 5, some brief remarks are made on 
the influence of body cross-sectional shape. 


NOTATION 
aspect ratio 
Cy wave drag coefficient, = Dwave/(qS) 
AC, drag-rise coefficient, =(D/qS)—(D/qS)u-o.s 
D total drag 
Dyave Wave drag 
M free stream Mach number 
q=4pV* 
S wing area 
slenderness parameter, span/length 
V free stream velocity 


A,,, angle of mid-chord sweepback 


p free stream density 


thickness/chord ratio 


2. Experimental Results 


Limitations of the sonic area rule can be investigated by analysing suitable 
experimental data. Possible sources of such data are wind tunnel and free flight 
tests. It is fortunate that some suitable test results do exist and these are 
discussed here. 


The wind tunnel data presented by Spreiter’’ and Page’ were obtained at 
Reynolds numbers, based on the mean wing aerodynamic chord, of one to two 
million. Page’s transonic tunnel results apply to delta wings and are more accurate 
than the transonic bump results, for rectangular wings, used by Spreiter. Thickness/ 
chord ratios varied from 2 per cent to 10 per cent and the aspect ratios from 4 to 6. 


The free flight results’*: *’ were obtained at a Reynolds number, based on the 
mean wing aerodynamic chord, of about six million. All body-alone drag results 
are considerably more accurate than the wing drag results, the maximum errors of 
which could be 25 per cent. The bodies had thickness/length ratios between 6 per 
cent and 8 per cent, while the wings of the wing-body combinations had 6 per cent to 
12:5 per cent thickness/chord ratios. Wing aspect ratio varied from 2 to 4. All the 
test vehicles were accelerated up to speed by 3 in. cordite rocket motors. During 
the subsequent coasting flight the total vehicle drag was derived by a double 
differentiation of the trajectory, as determined from three-station reflection Doppler. 
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Wing drag was determined from the difference between the drags of the wing-body 
combination and the basic body alone, interference effects on the drag 
being negligible. 


It must be noted that some accurate experimental drag results for swept-back, 
tapered wings are still desirable because of the rather poor accuracy of the free 
flight results'®’. 


3. Range of Applicability 


The range of applicability of the sonic area rule for rectangular wings has been 
investigated by Spreiter®’. If the drag-rise depends solely on the cross-sectional area 
distribution then, for a given family of wings, it can be shown'’ that [ACp/(Ar*)]y-, 
would be constant, independent of Ar'/*. Here AC, is the drag-rise coefficient 
based on the exposed wing area, A is the aspect ratio, - is the thickness/chord ratio 
and M denotes the free stream Mach number. A simpler derivation of this result 
can be obtained from the theoretical version’ of the sonic area rule. This version 
treats wave drag, not drag-rise. Since it is derived from supersonic theory it satisfies 
the first-order supersonic similarity law and so'’’ the wave drag parameter C,,/(A7*) 
must be a constant for a given wing family. Therefore, the sonic area rule would be 
expected to be applicable where or AC) /(Az*), is independent of Ar' 


3°5 


3:0 


fac\ 20 / 


1:5 rd 
/ / WINGS 
/ (REFERENCE 4) 


5 20 25 30 
Ar’? 
Fictre 1. Similarity analysis of (AC,/7°*),_, for wings with N.A.C.A. 63AOXX sections 
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Figure 2. Effect of transonic slenderness parameter and «7! * on the sonic wave drag 
parameter (AC,,/7°/*)y_.. 


For rectangular wings Spreiter®’ showed that (AC)/7°'*)w-, was a linear 
function of Az'* for Az'/* <1. Page‘ carried out wind tunnel tests on delta wings 
and obtained a wider range of linearity. These two sets of results are illustrated in 
Fig. 1. Clearly the scope of these results should be extended to include wings with 
swept-back trailing edges. Some drag-rise results for a series of swept-back, tapered 
wings have been obtained at W.R.E. using the free flight technique'*’; the range of 
linearity was found to become larger as the mid-chord shearback A tan A,,, 
increased and the wing became more “slender.” Taper ratio was not an important 
parameter. 


It is possible to obtain a much neater presentation of all this experimental data 
by introducing the slenderness parameter"? o (=span/length) and plotting 
against as in Fig. 2. Clearly (AC)/7*'*)w-, varies linearly with 
o7'!® (or Ar'!*) for all the wings provided that o7' * is less than about unity. This 
result is a generalisation of Spreiter’s result for rectangular wings and is likely to be 
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Figure 3. Similarity comparison of drag-rise results for wings and equivalent bodies. 


applicable to wings of arbitrary plan form. Hence it is suggested that the sonic area 
rule will be applicable to thin wing and slender body combinations when the wing 
slenderness parameter «7''* is less than unity. 


4. Application to Wings with Round-Nosed Sections 


Theoretical investigations of the drag of non-lifting wings with round-nosed 
sections show that a so-called leading edge drag force exists when the component of 
free stream Mach number normal to the edge is less than unity. This edge drag is 
incorporated in the wave drag, deduced from area rule considerations.* However, 
some additional leading edge drag would appear if the flow in the nose region were 
no longer of the reversible subsonic type, normal to the leading edge, but was 
essentially transonic in character. Then an equivalent body of revolution, having 
the same distribution of cross-sectional area as the wing, would have a lower 
transonic drag-rise than the wing. 


Whitcomb’s original statement of the transonic area rule''’ implies that the 


*The author is indebted to a referee for pointing out that Ref, 5 is incorrect in assuming 
that the edge drag should be added to this wave drag. The leading edge drag force is 
added only to that theoretical wave drag calculated from the integration over the wing of 
the product of local pressure and local surface slope. 
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drag-rise of a wing is approximately the same as that of its equivalent body. From 
the foregoing discussion it can be seen that this would be true only when the wave 
drag was very much greater than any additional drag-rise of the leading edge force. 
That there may exist configurations for which this is not true is shown by Fig. 3, 
which presents the drag-rises®’ of four wings and their equivalent bodies of 
revolution. These wings had round-nosed R.A.E. 104 sections, with Ar'* being 
0:841, 1:189, 1682 for Wings 1, 2, 3 (and 4) respectively. Also, the mid-chord 
shearback A tan A,,, was equal to 2 and the taper ratio was 0-408. The slenderness 
parameter o7''* varied from 0:50 to 0:96, being less than unity. Thus, from 
Section 3, the sonic area rule would be expected to be applicable. 


A systematic difference between wing and equivalent body drag-rise is evident 
in Fig. 3. The assumption that any flow separation effects (due to shock wave and 
boundary layer interaction) are unimportant in drag-rise comparisons has enabled 
estimates to be made of the additional leading edge drag force. This is found’ to 
be about twice the drag-rise of the theoretical leading edge force'*’. It is clear that 
flow separation effects may be incorporated in the estimate of the additional leading 
edge drag-rise. Unfortunately, these effects cannot be isolated, since no satisfactory 
method of allowing for them exists. 


Pending a more complete investigation of the effective leading edge drag force 
it is suggested that, near sonic speed, an additional leading edge drag-rise, equal to 
twice the theoretically predicted leading edge drag-rise, be incorporated in drag 
estimates of wings having thickness/chord ratios of about 6 per cent to 9 per cent. 


A wave drag estimate, obtained from the theoretical version of the sonic area 
rule’, is also shown in Fig. 3. It is not surprising that this estimate is very much 
greater than the sonic drag-rise of the equivalent body; the essentially supersonic 
theory of Ref. 2 would not be expected to give satisfactory wave drag estimates at 
sonic speed. 


5. Influence of Body Cross-Sectional Shape 


One group of free flight results’ referred to wing-body combinations having the 
same distribution of cross-sectional area, but different body shapes. The bodies had 


FiGurE 4. Typical test vehicle. (Maximum body waisting in the plane of the wing.) 
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thickness/length ratios of 8 per cent and so could be considered “slender.” A 
typical test vehicle, with stabilising fins, is shown in Fig. 4. The body shapes used 
were circular, elliptic with the major axis parallel to the wing plane and elliptic with 
the major axis normal to the wing plane. A constant total cross-sectional area was 
maintained in the region of the wing. 


The test results’’’ showed no significant variations in drag with body shape. 
Thus body cross-sectional shape has no significant effect on the transonic drag of 
these wing-body combinations. It is suggested that this result applies to all 
configurations with smooth bodies that are slender. For non-slender bodies, however, 
cross-sectional shape may be a significant variable. 
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The Stability of the Short-Period Motion of 
an Airframe Having Non-Linear Normal 
Force and Pitching Moment Curves 


P. A. T. CHRISTOPHER 


(Department of Aerodynamics, The College of Aeronautics, Cranfield) 


SuMMaRY: The analysis of the stability of second-order non-linear systems by 

Poincaré’s method of singular points in the phase plane is briefly presented 

and used to determine stability criteria for the short-period motion of an 

airframe having non-linear normal force and pitching moment curves. These 

results are then compared with those obtained by the application of quasi- 
linear stability theory. 


1. Introduction 


The linear theory of airframe dynamic stability is based on the assumption that, 
during a disturbance of the steady motion, powers and time derivatives higher than 
the first can be neglected (see Ref. 1). In some applications the linear theory is 
invalid, even as a first approximation, since the governing differential equations are 
markedly non-linear in form. An example of this kind arises in relation to the 
stability of low aspect ratio wings having a leading edge separation, the separated 
flow rolling up into approximately streamwise vortices which produce non-linear 
force and moment variations with incidence. Similar non-linear variations occur in 
wings at transonic speeds (see Ref. 2), on wings and bodies at hypersonic speeds and 
on most configurations at incidence near the stall. 


Except for those cases where it is possible to separate the variables or the 
equation is exact, no explicit general solutions to non-linear differential equations 
are known. Between the extremes of the exact solution and the approximate quasi- 
linear one a variety of higher approximations can be derived; these are usually based 
on analytic iteration techniques which yield series solutions of differing degrees of 
accuracy. For equations of the second order and lower the “phase plane” method, 
which is qualitative, will give considerable insight into the nature and stability of 
solutions and can, by the use of Liénard’s or other constructions, give quantitative 
solutions. In Ref. 3 the phase plane method is discussed and applied to examples 
taken from elementary mechanics, structural engineering and electrical machine 
theory. 
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The longitudinal free motion of an airframe is often expressed in terms of a 
differential equation of at least the fourth order; the principal modes of oscillation, 
the “short peroid” and “phugoid”, are then interrelated (see Ref. 1). In some 
applications, for example many guided missiles, the short-period motion, which is 
governed by an equation of the second order, is of dominant interest. When the 
influence of non-linear normal force and pitching moment variations with incidence 
are included in this equation, it provides an interesting example of the phase 
plane method. 


NOTATION 


a,b,c,d constants in Poincaré’s theory of singular points 
f,.f. functions defined in equation (17) 
m airframe mass 
q angular velocity about axis of pitch 
time 
v=dx/dt 
w velocity along axis of yaw 


x any independent variable, see equation (1) 


aM) 


on n=0 

max (2M) 
/ 


| constants in the force and moment relations of equation (12) 


A, =(U,+ Me — 
A,=(U,+2,) m,—m,2Z; 

B moment of inertia about the axis of pitch 
B,=(U, + 2.) Ma + mM, + 
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B,= 32, 
D the operator d/dt 
G=-2A, 


B, A, 
B, Ay 


moment about axis of pitch 


H=B,(1- 


forcing function of the elevator 

velocity along the axis of roll (longitudinal body axis) 
force along axis of yaw 

elevator angle 

angle of elevation 

root of characteristic equation (9) 


ordinate of w referred to singular point 
2. The Phase Plane Method 
The method is concerned with equations of the form 


d*x 
dt’ 


or degenerate forms. Writing v=dx/dt, then d*x/dt®? =vdv/dx and (1) becomes 


dx 
+ B(x) +C (x)=0, 


dv 
v +B(x)v+C (x)=0. 


The graphs of the solutions of this equation in the xv-plane, known as the phase 
plane, are referred to as “integral curves” and through each ordinary point in the 
plane there passes only one such curve. 


Equation (2) may be rewritten as an equivalent pair of equations 


dv 
dt = —C(x)—-B(x)v 


which define a field vector having components dv/dt and dx/ dt; this vector is always 
tangential to the integral curve and indicates the direction in which f is increasing. 
When B (x)=0, and C (x) is an anti-symmetric function, the integral curves are closed 
and the slope of the field vector is 


_ dv —C(x) 


2 
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The variables can be separated and, provided that the integral is soluble, the integral 
curves may normally be expressed in analytic form. The resulting motion is 
periodic, this being characteristic of conservative systems. For the general case, 
explicit integration is not possible and solutions may be obtained by geometric 
constructions, such as that of Liénard. To examine the stability of the solutions of 
(2) or (3) it is not necessary to have these explicitly; for this purpose use may be 
made of Poincaré’s theory of singular points, which is now briefly summarised. 


3. Singular Points 


The stationary positions of equilibrium of equations (2) or (3) correspond with 
the singularities of the equivalent equation 


dx v 
and analysis of the character of these singularities gives considerable insight into the 
nature of the motion near these points. More generally, consider the singularities of 
the equation 


dv P(x,v) 
dx (x, v)” (6) 
which are defined by P(x, v)=Q(x,v)=0. Since the origin can always be changed 
to correspond with the singular point, then analysis can be restricted to singularities 
at the origin. 


When (6) has a singularity at the origin then it is assumed (Poincaré) that it may 
be written in series form 


dv ax+bv+p(x,v) 
dx cx+dv+q(x,v) 
where p(x, v) and q(x, v) are the remaining terms of series whose lowest terms are 
of second degree at least. Further, if the constants obey the inequality 


A= =ad—bc #0, 


la 

c d 
then the integral curves behave, in the neighbourhood of the singularity, as if p (x, v) 
and q (x, v) were absent. 


The singularities of the reduced equation 


dv ax+bv 


are of four distinct types. 
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Ficure 1. Types of singularities. 


The four types of singularities are 


(i) nodal points (or nodes), which are either stable or unstable, (a=d=0, 
c=1,6>0 or d=0, a=b=c=1) 


(ii) centres, which indicate periodic motion, (b=c=0, d=1) 
(iii) spiral points, which are either stable or unstable, (a=1, d= —1, b=c #0) 
(iv) saddle points, which are unstable (a=d=0, c=1,b <0). 


The nature of the integral curves in the region of these singularities is shown in 
Fig. 1. If the singularity does not correspond with one of the four cases (i)-{iv), then 
it may be reduced to one of these forms in a new variable by means of non-singular, 
linear transformations. 
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Equation (8) may be replaced by the equivalent pair of equations for the 
components of the field vector 


=t?=ax+h 

8 
(8a) 
at =x=cx+dv, 


Also, x and v may be transformed into x, and v, by the linear transformation 


x,=ax+ Bu 


v,=yx+ dv, 


implying that the integral curves are now referred to new oblique co-ordinate axes. 
A further transformation is now sought which makes the slope of the field vector in 
the x,v,-plane 


dv, _ 
dx, = rx, ° (9) 


This form is directly comparable with certain of those in (i)-(iv) and it is left to 
obtain A, and A, in terms of the coefficients in equation (8). Since the transformation 
is to be non-singular, then v, 40, x, 0 for all values of x and v, and thus 
(25—y8)~0. Equation (9) may be written in the equivalent form 


dt =0,=A,2, 
(9a) 
dx, 
"at =X, =A,xX, 
and, since and 7,=yx+ 4, then, upon substitution for x, %, x,, 
x, and 7%, in (9a), there arise the equations 
a (cx + dv)+ B (ax+ bv)=A, (ax+ Bv) (9b) 
y (cx+ dv) + 3(ax+ buv)=A, (yx+ bv) 


Since these must be true for all values of x and v, then the coefficients of x and v 
must vanish (i.e. the coefficient may be equated) giving the equations 


a(A, —c)—Ba =0 \ 


(10a) 
—ad +B(A,-d)=0 | 


and y(A,-—c)— 8a 0 (106) 


~yd + 8(A,—b)=0. 
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The coefficients z, 2, y and 6 are not all zero; therefore (10a) and (105) have a 
common eliminant 


| (A-c) 
-d  (A-b) | 


or A?—A(b+c)—(ad—bc)=0. . (0c) 


Criteria for distinguishing the type of singularity are obtained from the “character- 
istic equation” (10c) which has the roots 


A,.3=4 {(b+c)+[(b +c)? + 4 (ad — bc)}'/*} 


and, since ad—bc 0, then A can never be zero. Three important cases are 
distinguished, as follows. 


Roots real and unequal. These arise when the discriminant is positive. The 
sub-cases are: if A > 0 then A,/A, is negative and corresponds to a saddle point, and 
if 4 <0 then A,/A, is positive and corresponds to a nodal point. In the second 
case, if A, and A, are both negative the node is stable, whereas if A, and A, are both 
positive the node is unstable. 


Roots complex conjugate. Complex roots occur when the discriminant is 
negative. If the real part of the root is negative the integral curves are stable spirals, 
whereas if the real part is positive the curves are unstable spirals. When the roots 
are purely imaginary the singularity is a centre; however, under these conditions the 
singularity of (8) is not necessarily that of (7) and the higher order terms in p (x, ”) 
and q (x, v) have to be considered. 


Roots real and equal. The singularity is a node and its stability is governed in 
the same way as if the roots were unequal. 


4. Equations of Longitudinal Motion of an Airframe 


In analysing the longitudinal dynamic stability of guided missiles it is customary 
to consider only the short-period motion, the velocity of flight being assumed 
constant during a small disturbance. In this case the linear equations of motion 
become (see, for example, Ref. 5) 


(D-Z«) w—-(U, +2) = 2 | 
—(m,D+m,)w+(D—m,) 


(11) 


For the present problem it is intended to include non-linear variations of the 
normal force Z and pitching moment M with vertical velocity w and to introduce 
these effects into equations (11). It is assumed that the force and moment contribu- 
tions due to w, written Z (w) and M (w) respectively, may be expressed in any of the 
graphical forms shown in Fig. 2. The important quality of these curves is that they 
are anti-symmetrical in form, characteristic of most symmetrical configurations. The 
anti-symmetric character allows the forces and moments to be written as power 
series in odd powers of w. 
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-Z (w) m(w) 
I¢ 
\ 
4 
4 / \b 
a 
Figure 2. Normal force and pitching moment curves. 
Thus 2) +... 
(12) 
M (w) 
B w+... 


When these non-linear expressions are substituted into the equations of motion, 
they become 


. .)+(D—m,) . (14) 


Eliminating 6 between (13) and (14) then gives 


(D—m,) (Ww —z.w-z,w'—. . mw + +...) 


=[z,(D—m,)—m,, (U, + 2,)} 


or, upon collecting terms, gives 


—[(U, + 2,) ma + my + Zw + W—[(U, + Me Zw] W- 


. .=[z,D—-(U,+2,) m,—2,m,] (15) 
This equation may be written in the form 
w—f,(w)w—f, (w)=Q (n) (16) 
which, for the case of the free motion, Q (n)=0, reduces to 


and is of the same form as (1). 
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Note: It is important to appreciate that in a non-linear system the principle of 
linear superposition no longer applies, and the solution of (17) is no longer a 
complementary function to the complete solution of (16). In linear theory, if the 
corresponding equation to (17) can be shown to be stable, then the corresponding 
equation to (16) will also be stable provided that Q (7) is of finite magnitude. The 
stability of (16) and (17) are not related in this simple manner. 


The variable w cannot be eliminated between equations (13) and (14) in such a 
way as to give an equation in 4 which is of the same form as (16) and therefore the 
stability of the 6 motion cannot be obtained by the method of singular points. 


5. Stability of the Short-Period Motion 
Writing w =v, then (17) becomes 


dv f.(w)+f, (w)v 


dw v 


‘ ‘ (18) 


which has singularities when f, (w)=0, 7 =0. The number of singular points depends 
on the number of unequal roots of the equation 


f.(w)=[((U,, + 2,) mo — w+ [((U, + 2.) — +... . 
[(U + 2q) Mon—, — M Zon] +... .=0, (19) 


and therefore on the number of terms, m, used to represent the force and moment 
curves. For simplicity only two terms are considered in the remainder of the 
analysis, although this can readily be extended to any reasonable number of singular 
points; hence 


f. (w)=[(U, + 2,) Me —M Ze] ((U, + Z,) —m,z,] =0. (20) 


When the coefficients of w and w* in (20) are the same sign, then only one singularity 
exists and this is at the origin, while if the coefficients are of different sign singular 
points occur at 


w= + { —(U,, + Me + M Ze 1/2 


v 
(U, + Z,) ms — v=0, . 


as well as at the origin. 
Consider the singularity at the origin and write (18) in the first approximation as 
dv A,w+B,v 


dw 


where A, =(U,+ 2%) Me — 
B,=(U, + 2%) My + + Zu, 
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terms of w* and higher degree being neglected. Comparing with the standard form, 
equation (8), a=A,, b=B,, c=0 and d=1, and the roots of the characteristic 
equation are 


{(B, + (B,?+4A,)'?}. 
The different cases are :— 


B, <0, A, <0 such that B,?+4A,<0. The roots are complex conjugate 
and the singularity is a stable spiral point. 


B, <0, A, <0 such that B,?+4A,>0. Both roots are real and negative, 
corresponding to a stable node. 


B,<0,A,>0. The roots are of opposite sign, giving rise to a saddle point. 


From the definition of A, and B,, it can be seen that the sign of A, is largely 
determined by that of m,, and that of B, by the sign of z,. In normal circumstances 
Z» Will be negative and so will B,. As a first approximation m,, will be proportional 
to the centre of gravity margin and thus A, may be interpreted as a static stability 
parameter. A, <0 implies a large amount of static stability which, when combined 
with B, < 0, the latter condition implying positive damping in the linear case, gives 
a damped oscillation. Increasing A, corresponds with decreasing the static stability, 
A, =0 being the boundary between a stable node and saddle. 


The equation of free motion for the linear case is obtained from (15) by 
neglecting second and higher powers of w and the terms in », giving 


. . . 


For stability of the linear equation both B, and A, must be negative, which is the 
same criterion as that for a stable singularity at the origin in the case of the 
non-linear equation. 


The singularities away from the origin may have their co-ordinates written as 
w=w,, v=0 and, referred to these co-ordinates as origin, w may be written 


w=Weté . ‘ ‘ (24) 


Equation (18) now becomes 


dé Vv (25) 
which, upon substitution for f, and f,, gives 


dv A, (w,+ A, (wy + + [B, + B, 


dé v 


where A, =(U,+z,) m,—m,z, and B, = 3z,. 
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Expanding, 


dv _ [A.w.+ Awe 3AweE+. +B, (wet 
v 
Now A,w,+A,w,*=f, (w,)=0, 


and therefore in the first approximation 


v 


dv _ G£+Hv 


ad dé v 


where G=A,+34,w,?=A,+34, (- “) 
3 


| A 1 B, A i 
and H=B,+B,w,2=B,-B, (2) =B,|1- 
In order to produce this singularity, A, and A, must be of opposite sign. Comparison 
with the standard form, equation (8), in which x is replaced by &€, gives 
a=G, b=H, c=0 and d= 1 and the roots of the characteristic equation are 


A,,2=4 {H+(H’? + 4G)'?} 
The different cases are : — 


H <0, A, <0,A,>0. Here G is positive and the singularity is a saddle 
point. 


H <0, A, >0, A, <0. Now G is negative. When H*?+4G <0 the 
singularity is a stable spiral point. When H*?+4G => 0 the singularity is a 
stable node. 

The damping of the motion near the origin is dependent on B,, which in normal 
circumstances is negative, whereas in the present case H can change sign. Since 
A,/A, is always negative if the singularity is to exist, then B, must be positive if H 
is to change sign. 

Two further cases emerge : — 


H>0,A,<0,A,>0. Here G is positive and the singularity is a saddle 
point. 

H>0,A,>0,A,<0. G is negative. When H*?+4G <0 the roots of 
the characteristic equation are complex with positive real parts and 
correspond to an unstable spiral in the phase plane. When H*+4G >0 
the roots are real and of different sign and correspond to a saddle point. 
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dv _ (A, +3A,w,’) £+(B, + Byw,’) v 
dg 


FINN 
SN 
(b) A,>0,a,<0,H<0, 
Cy 
S Zz 
WY, (c)A,>0,A,<0,H>0. 
WY H*+4G<0 
© 
ZF 
Ficure 3. Integral curves in the region of the singularities, B, < 0. 
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It is now appropriate to plot the integral curves in the phase plane. In the 
cases where A, and A, are of the same sign, then the singularity occurs only at the 
origin and the non-linear terms do not change its character. When A, and A, are 
of different sign, then two additional symmetrically disposed singularities arise 
whose natures depend on the signs of H, A,, A, and the discriminant H?+4G. The 
three interesting cases are those where the singularities away from the origin differ 
from those at the origin and are sketched in Fig. 3. 


The curves of Fig. 3(a) correspond to an M (w) curve such as c in Fig. 2, and 
imply that for initial values of w between +, the airframe will perform a damped 
oscillation. When z,, and hence B,, is large, then B,7+4A,> 0 and the motion 
near the origin degenerates into a subsidence. In the quasi-linear approximation 
m, will change sign at the minimum in the M(w) curve, and (17) will become 
unstable when df, (w)/dw changes sign. 


Differentiation gives 


=(U, + 2) My — + 3 {(U, + 2) m, —m,2;} w’, 


which changes sign when 


Since the dominant terms in (21) and (27) are those containing m,. and m,, then it 
follows that the range of w over which damped oscillations will result from an initial 
disturbance is greater by a factor of approximately 3 than that predicted by the 
quasi-linear approximation. Any changes of sign of H due to the influence of 
B,=3z, have no influence on the nature of the singularities at + w,. 


The curves in Figs. 3(b) and 3(c) correspond with an M (w) curve such as 6 in 
Fig. 2. When H < 0, which it normally will be unless B, is large and positive, such 
that H?+4G <0, then any initial disturbance will cause the airframe to perform 
damped oscillations about +w,. From the quasi-linear theory df, (w)/dw changes 
sign at a value of w given by (27) and stability is indicated for larger values of w, 
provided that B, is negative. If B, is positive this implies that f,(w) will also 
change sign (this corresponds to the stalling of the airframe), which occurs when 
f, (w)=0 or 


1/2 
32; 
If the value of w given by (27) is smaller than that given by (28), there will exist a 
region between these values for which stability is predicted by quasi-linear theory, 
whereas if the converse is true no stable region exists. The non-linear theory 
indicates that if B, is such that H < 0, then the motion will always be stable in the 
sense that it performs a damped oscillation about +w,, whereas the quasi-linear 
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theory indicates stability about any value of w greater than that given by (27). If B, 
is large and positive such that H > 0, then non-linear theory indicates instability for 
all values of w, while quasi-linear theory indicates a stable region for values of w 
lying between w,, and w., provided that w,, < w,. 


6. Conclusions 


The important conclusions which may be drawn from this analysis are as 
follows. 


(1) The method is only suitable for a partial determination of the stability of 
the short-period motion of an airframe with the elevators undeflected. It is partial 
in the sense that it can only consider the influence of the non-linearities on the 
equation governing the vertical velocity w (or incidence) and not that involving the 
angular rate of pitch. In the non-linear problem these two equations differ in form, 
whereas in the comparable linear problem they do not. 


(2) Comparisons with quasi-linear theory (i.e. linear theory based on the local 
slope of the Z (w) and M (w) curves) indicate the following : — 


(a) The criteria for stability at low incidence is the same for both methods. 


(6) With airframes having pitching moment curves such as b and c in Fig. 2, 
the stability criteria and motions predicted by the two methods differ 
considerably at moderate to high incidence. With curve c, quasi-linear 
theory predicts instability when w has values outside the range between 
the maximum and the minimum, whereas non-linear theory indicates a 
stable region which is greater by a factor of approximately /3. With 
curve b and a Z(w) curve whose gradient does not decrease, quasi-linear 
theory predicts stability when w has values outside the range between the 
maximum and the minimum; whereas non-linear theory indicates that the 
only possible motion of the airframe is a damped oscillation about values 
of w=w, given by equation (21). With curve b and a Z (w) curve having a 
decreasing slope, which includes the case of the stall, quasi-linear theory 
indicates a stable region for values of w between those corresponding to 
zero stiffness w,, and zero damping w.; non-linear theory indicates 
instability for all values of incidence. 
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Summary: A set of design charts is presented for the calculation of transient 
temperature and thermal stress distributions in thermally thick plates subjected 
to aerodynamic heating. 


The method is particularly useful for determining temperatures and thermal 
stresses in plates with an arbitrary variation of the heat transfer coefficient and 
the adiabatic wall temperature of the boundary layer. The present method is 
based on repetitive applications of the exact analytical solution to a unit 
triangular variation of the adiabatic wall temperature and a constant heat 
transfer coefficient. The actual variation of the adiabatic wall temperature is 
represented as a series of straight lines while the heat transfer coefficient is 
approximated by a step function. The temperature distribution through the 
plate is separated into linear and “self-equilibrating” temperature distributions 
to facilitate thermal stress calculations; these distributions can be obtained 
directly from the design charts presented in this paper. 


The general principle of this semi-numerical method is also applied to 
thermally thin plates subjected to arbitrary heating conditions. 


1. Introduction 


In heat transfer calculations for plates subjected to aerodynamic heating or 
cooling it is often necessary to consider temperature gradients through the thickness 
of the plate. Whether these gradients are important depends on the magnitude of 
the heat transfer coefficient h at the exposed surface, the plate thickness d, and the 
thermal conductivity k. These three variables can be combined into a single non- 
dimensional parameter usually described as the Biot number B=hd/k. It will be 
assumed here that for B >0-1 the heated (or cooled) plate can be defined as 
“thermally thick”, which is a plate in which temperature gradients through the 
thickness are not negligible, while for B < 0-1 the plate can be defined as “ thermally 
thin”, in which temperature gradients are negligible. It should, however, be 
emphasised that this division is somewhat arbitrary and it is only introduced to 
facilitate practical calculations. 
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The one-dimensional heat flow by conduction in a medium whose thermal 
properties do not vary with temperature is governed by the partial differential 
equation 


eT 1a 


where T is the temperature, x is the distance, « is the thermal diffusivity and ¢ is the 
time. In problems involving aerodynamic heating of plates the boundary conditions 
are usually taken as 


and oT x=0, ‘ ; . (3) 


where k is the thermal conductivity, A(t) is the heat transfer coefficient and T,, is 
the adiabatic wall temperature. In the first boundary condition the rate of heat 
absorption by the material at the heated surface at x=d is equated to the rate of 
heat’ input by forced convection (aerodynamic heating), while the second condition 
assumes that the surface at x=0 is insulated (zero heat transfer). 


The heat transfer coefficient h is usually evaluated as a function of time, 
depending on the variation of aerodynamic conditions just outside the boundary 
layer. For more accurate calculations, however, the variation of 4 with the surface 
temperature should also be included, particularly in heat transfer calculations for 
turbulent flow”. 


The arbitrary variation of h(t) and 7, (t) in practical design problems rules out 
any possibility of a successful analytical solution for the temperature distribution; 
only solutions to relatively simple cases can be obtained. Solutions to the problems 
with a constant heat transfer coefficient and a suddenly applied temperature 7,,, are 
known and furthermore the case of arbitrary variation of T,, can be solved using 
Duhamel’s integral. Kaye and Yeh obtained a solution for the temperature 
distribution assuming that the heat flux through one surface varies as a quadratic 
function of time. Chen represented the heat transfer coefficient and the adiabatic 
wall temperature by constants and exponential functions and solved the differential 
equation by means of the Laplace transformation method. 


Several objections may be made to the use of the analytical methods such as 
those given in Refs. 2 and 3. Firstly, solutions may not be sufficiently realistic, due 
to the various simplifications which are necessary to get any solution at all; 
secondly, the analytical expressions are usually extremely complicated for practical 
applications; and thirdly, the analytical solutions usually converge poorly near t=0 
unless special asymptotic solutions are employed for small values of time. 


The numerical methods are based essentially on step-by-step integration of the 
heat conduction equation. These methods have the additional advantage over 
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analytical methods that the variation of thermal properties of the material with 
temperature and the temperature dependence of the heat transfer coefficients can be 
included. Similarly any radiation losses and other non-linear effects, e.g. free 
convection, can be included in the solution. Numerical methods applied to heat 
transfer problems are, in general, time-consuming and require relatively small 
intervals of time. On the other hand, they can be programmed for automatic 
computation and several such programmes for digital computers are already 
in existence. 


The semi-numerical method described in this paper is based on repetitive 
application of the exact analytical solution to the response to a triangular variation 
of 7, over a relatively large time interval. This is possible because of the linear 
character of equations (1), (2) and (3). The heat transfer coefficient h is assumed 
here to be constant in any given time interval, while the variation of T., is repre- 
sented by a series of straight lines. A similar approach has been used in Ref. 4, 
where the adiabatic wall temperature 7,,. was also represented by a finite series of 
unit triangles. The main advantage of this method is that the computing time is 
appreciably reduced, compared with standard numerical methods, if special charts 
are available for the temperature distribution due to linear variation of T,,. Such 
charts are presented in this paper and they can be used for calculating either the 
temperature or the thermal stress distributions. 


To facilitate calculations involving thermal stresses, the temperature distribution 
through the thickness of a plate is separated into linear and “self-equilibrating” 
temperature distributions”; the linear distribution contributes to the overall 
deformation, either curvature or expansion in the plane of the plate, while the self- 
equilibrating distribution produces zero overall deformations. 


NOTATION 


T (x, t), T (z,t) temperature distribution 
T(n,N) temperature distribution (in non-dimensional distance and time) 
T.. adiabatic wall temperature 
T, adiabatic wall temperature at time N, 
Ta (N) mean temperature 
AT (t), AT (N) temperature gradient 
T*(z,t),T*(n,N)  self-equilibrating temperature distribution 

x distance measured from the insulated surface 
z=x-—d/2, distance measured from the middle surface 
n=x/d, non-dimensional distance 
d_ plate thickness 
t time 
x=k/(pc), thermal diffusivity 
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k thermal conductivity 
h_ heat transfer coefficient 
N=xt/d*, non-dimensional time 
6 non-dimensional time interval 
N,=r6 
r,j integers 
B=hd/k, Biot number 
= ad? 
a rate of increase of T,,, defined by equation (11) 
E Young’s modulus 
a coefficient of thermal expansion 
Poisson’s ratio 
@ thermal stress 
p density 
c specific heat 
w, j positive root of w;tanw;=B 
F(N,, B,) temperature functions for B= B, and N=N, 
B, Biot number at N, 
R(N,) temperature response at N, 


T =\ Te~?* dN, transformed temperature (Laplace transform) 


p function in the Laplace transform 
exp(_) exponential function 
erfc( ) complementary error function 
i’ erfc( ) repeated integral of the complementary error function 
erfc (0)=2-'/U (1+r/2) 


) gamma function 


2. Temperature Distribution, Self-Equilibrating Temperature 


The solution of equation (1) can be separated into linear and self-equilibrating 
temperature distributions, i.e. 


T (2, 1)=Tm(t)+AT +T* (2,0), 4) 


is the distance measured from the middle surface towards the heated surface at x=d, 
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T.» (t) is the mean temperature, AT (t) is the linear temperature gradient and T* (<, t) 
is the self-equilibrating temperature distribution, for which 


a/2 
-a/2 -d/2 


Thus only the first two terms on the right hand side of equation (4) contribute to the 
overall deformation of the plate. 7,, and AT terms produce thermal stresses in the 
presence of edge restraints. On the other hand, stresses due to T* occur even if all 
edges are completely free. 


Using equations (4) and (6), it follows that 


| Taz 


and AT=-, « «» « & 


Alternatively, AT could be obtained by equating the first moment of the temperature 
T (z, t) to the first moment of the AT (1) z/d distribution about the middle plane. 


For a plate with all edges free, the thermal stress distribution through the 
thickness of the plate, away from the edges, is given by 


/ 


d/2 d/2 


/ -d/2 


Obviously, the stresses given by equation (9) fall to zero in a narrow “boundary 
layer” along the edges. 


Now, using equations (4), (7) and (8), it follows that the thermal stresses given 
by equation (9) are proportional to the self-equilibrating temperature distribution 
T*, i.e. 


The mean temperature 7,, induces thermal stresses only if one or more 
edges are restrained against expansion in the plane of the plate. In such cases 
compressive or tensile stresses are induced if the plate is respectively heated or 
cooled more rapidly than its surrounding structure. The effect of the temperature 
gradient AT is to produce a uniform curvature in the plate and therefore, in the 
absence of any edge restraints, the plate bends into a spherical surface. Because of 
the restraint along the edges, which prevents the edges from leaving the plane of the 
supports, the temperature gradient AT produces certain reactions along the 
boundaries and bending stresses in the plate. Ref. 1 gives a full discussion of 
thermal stresses due to temperature components T,,, AT and T*. 
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3. Linearly Varying Adiabatic Wall Temperature 7,., 


If the adiabatic wall temperature varies linearly with time and the initial 
temperature is assumed to be zero, then 


where a is the rate of temperature increase. The assumption 7,,,=0 at t=0 does 
not affect the results because it is always possible to measure 7,,, and hence the 
temperature distribution in the plate, from some initial equilibrium temperature. It 
is shown in Appendix I that, for a constant heat transfer coefficient A and a linearly 
varying 7,,,, the temperature distribution is given by 


T(n,N) 


Sin w; COS 


7 _w2Nn 


where the w, are the positive roots of the transcendental equation 


The parameters ?, n, N and B used in equation (12) are defined as follows : — 


N=kxt/d*? (non-dimensional time) ; ; (16) 
B=hd/k (Biot number) ; 


Using now equations (4), (7), (8) and (12), it is found that 


-12(n-9 (= exp(—w2N). (20) 


In general, the complete temperature distribution through the thickness of the 
plate is not required and for design purposes only the following quantities are of 
practical interest: surface temperatures T(0,N) and T(1,N), mean temperature 
T,.(N), temperature gradient AT(N) and the self-equilibrating temperatures 
T* (0,N) and T*(1,N). These component temperatures have been plotted in 
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4G Next/d" 


12 20 24 
Ficure 1. Surface temperature at n=O (insulated surface). 
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Ficure 2. Surface temperature at n=1 (heated surface). 
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0 0-4 08 20 24 28 
Ficure 3. Mean temperature T,, (N). 
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Ficure 4. Temperature gradient AT (N). 
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Ficure 5. Self-equilibrating temperature at n=0 (insulated surface). 
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Ficure 6. Self-equilibrating temperature at n=1 (heated surface). 
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Figs. 1 to 6 as functions of the non-dimensional time N = «t/d? for a range of values 
of the Biot number B=hd/k. In general, for B < 0-1 the plate can be regarded as 
thermally thin, for which any temperature gradients through the thickness are 
negligible. For such cases, only the mean temperature T,, is required and it can be 
calculated by means of equation (26) given in Section 5. 


For small values of the non-dimensional time N, asymptotic expressions derived 
in Appendix II can be used for various temperature components. 


4. Arbitrary Variation of T,,, 


In practice T,,, can be represented by a series of straight lines constructed by 
contributions from elementary triangular variations, as shown in Fig. 7. It is easily 
seen that the triangular variation is obtained from three separate linear variations 
starting at times N equal to N,_,, N, and N,,, respectively. This triangular variation 
may be regarded also as a “temperature impulse” applied to the heated surface of 
the plate. The magnitude of 7... at the time N, is denoted by 7, and the average 
heat transfer coefficient in the time interval N,_, to N,,, is assumed to be given by 
h,. Thus the corresponding value of the Biot number can be calculated from 
B,=h,d/k. For simplicity it will be assumed that N, —N,_, =4 (constant), or N, =r; 
but it should be emphasised that unequal intervals of time could also be used. 


| 


T, 
T, A 
| \ \ 
T, | \ \ \ 
for T,,, by a series of \ - | | { | { | 
triangular variations. | : = Noy 
| | | \ \ 
° N, Ng Ng Ne Nees 


Assuming that any of the non-dimensional temperature parameters represented 
on the design charts in Figs. | to 6, i.e. T (0, N)/®, T (1, N)/®, etc., is denoted by 
the function F (N, B), which is a known function of N and B, then the temperature 
response R (N) to the application of a linearly varying T,, after the time N, is given 
by (see Figs. 7 and 8) 
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Ficure 8. Triangular variation of adiabatic wall temperature. 


After the time N, has elapsed 


R(N,)= F(N,, B,)- = F(N,, B,)+ B). . 
while for N, 
R(N,)= F(N,, B,)- F(N,, B,)+ F(N,, B,)+ F(N,, 
- 21 F(N,, B,)+ F(N,,B,). . (23) 
For the time N, the general formula becomes 
8R(N,)=T, [F(N,, B,)—2F (N,_,, B,) + F + 
+T,(F(N,_,. B,)—2F (N,_., B.)+ F B,)] + 
+T, [F (N,_., B,)—2F (N,_3. B,)+ F + 
+ T,_. (F (N;, (N:. B,.)+ F (N,. B,-.)] + 
+T,_, B,,)—2F (N,, B,-,)] + 
+T,[(F C(N,, B,)]. (24) 
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Fquation (24) can be arranged in a tabular form and any required temperature 
components can then be computed rapidly with the aid of the design charts shown 
in Figs. 1-6. The method is sufficiently accurate for engineering purposes, even if 
large time intervals are used. 


It should be emphasised that 7,, T,,. . . T, can be either positive or negative, 
depending on whether the case considered represents acceleration or deceleration. 


If the adiabatic wall temperature 7,,, increases linearly from 7,,=0 to 7, in 
time N,, thereafter remaining constant, then this variation can be represented as a 
difference of two linear variations, as shown in Fig. 9. If, in addition, the Biot 
number B does not vary with time, then the response function can be simply 
expressed as 


and R(N)= [F(N, B)-F(N-N,,B)); N>N,. 


5. Thermally Thin Plates 


It is usually assumed that for B < 0-1 plates can be regarded as thermally thin 
and any temperature gradients across the thickness are negligible. It can be shown 
that for linearly varying adiabatic wall temperature T,,=at, the mean temperature 
T. of a thermally thin plate is given by 


Tm (N) 
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For large values of N equation (26) reduces to a straight line variation 


apcd 


or h 


(27a) 


where apcd/h represents the temperature difference between the adiabatic wall 
temperature and the mean temperature in the steady state. 


The mean temperature in a thick plate for large values of N can be obtained 
from equation (18). 


Hence 
aped,,_ , 
or h (1+4B). . (28a) 


It is therefore clear that equation (26) overestimates the mean temperature, if used 
for thick plates. Equations (27a) and (28a) become identical only for B=0. 


Equations (26) can be used for calculating the mean temperature due to a 
triangular variation of 7, and hence the semi-numerical method described in 
previous sections is also applicable to thin plates. 
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Appendix I 


SOLUTION OF THE DIFFERENTIAL EQUATION OF HEAT FLOW 


The basic equation (1), together with the boundary conditions (2) and (3) and linear 
variation of the adiabatic wall temperature, can be expressed in non-dimensional form as 


a(T/®) 


“an? ON 
0(T/®) _ BIN-(T/)) ; n=i (30) 
On 
On 
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and T/P=0 at N=0, (32) 


where the last equation implies zero temperature at N=0. In practical applications, 
however, the initial temperature at N=0 is constant and therefore the temperature 
distribution T (n, N) represents temperature increase above the initial temperature. 


Introducing the single-tailed Laplace transform 


~ 


equations (29), (30) and (31) can be transformed into 


@(T/®) _ 


dn? (4) 
d(T/®) 
dn 
and 
dn 


where equation (32) has been used in deriving equation (34). The solution of equation 
(34) with the boundary conditions (35) and (36) is given by 


cosh p'/?n 
T/D)= - 37 
(7/9) [((p'/?/ B) sinh p'/?+cosh p'!/?] G) 
Applying now the inversion theorem to equation (37) it can be shown that 
T (n, N) te sin w, COS Wyn 
2> —wZN). . (12 
+ 4(n? —1—2B-')+ (w, bain w, 008 w,?N). (12) 
where the w, are the positive roots of the transcendental equation 


The first six roots w, of equation (13) are given in Ref. 5 for a range of values of the Biot 
number B. 


Appendix II 
SOLUTION FOR SMALL VALUES OF NON-DIMENSIONAL TIME 


For small values of the non-dimensional time N, the parameter p in equation (37) 


must be large“ and hence sinh p'/? ~ cosh p'/? ~ 4expp'’*. The values of T/® at 
n=0 and 1 from equations (37) can therefore be approximated by the following 
equations. 
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T  2exp(-p') 

@ ~ p*(1+p' 2B-1) ’ n=0 ‘ (38) 

and ~ p? n=l. ‘ ‘ (39) 


In view of the invariance properties of the Laplace transform with respect to the second 
variable (distance), equations (7) and (8) can be used for calculating the transformed mean 


temperature T,, and temperature gradient AT. Thus, substituting equation (37) into (7) 
and (8), it follows that 


1 


Using now the table of Laplace transforms from Ref. 5, it can be shown that 
2B-? { Brany” it erfc (sy) 
—exp (B+ B?N) erfc (5 (42) 
=B-? { ly (4NY'? erfc (0)—exp (B?N) erfc (BN'/*) \ (43) 
=B-* { (- Br (4NY’? i” erfc (0)+ exp (B°N) erfc (BN'!?) (44) 
AT) +6 (B 2. (—1y~-' (4NY’? erfe (0)+ 
+exp (B°N) erfe (BN!) } +24B-* {5 (- 1yB’ (4Ny'? rerfe 
—exp (B+B?N) erfe (sy +BN*!?) } 


where i’ erfc ( ) is the r repeated integral of the complementary error function erfe ( ). 


The self-equilibrating temperature components can be calculated now from 
equations (42) to (45) and equation (4). Hence 


T* (0, N)/P={(T (0, N)—T,, (N)+44T (N))} /® (46) 
and 7*(1,N)/®=({(T (1, N)—T,, (N)—44T (N))} /®. (47) 
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No definite rule can be laid down for the maximum value of N up to which the 


approximate expressions can be used instead of the exact expressions given by equations 
(12) to (20). If the approximate expressions are used they should always be checked 
against the exact values for the maximum value of the non-dimensional time N likely to 
be used in the calculations. 


REFERENCES 


PRZEMIENIECKI, J. S. Design of Transparencies. Journal of the Royal Aeronautical 
Society, Vol. 63, pp. 620-636, November 1959. 


Kaye, J. and Yew, V. C. M. Design Charts for Transient Temperature Distributions 
Resulting from Aerodynamic Heating at Supersonic Speeds. Journal of the Aeronautical 
Sciences, Vol. 22, pp. 755-762, 786, November 1955. 


Cuen, S. Y. Transient Temperature and Thermal Stresses in Skin of Hypersonic Vehicle 
with Variable Boundary Conditions. Trans. A.S.M.E., Vol. 80, No. 7, pp. 1389-1394, 
October 1958. 


Hit, P. R. A Method of Computing the Transient Temperature of Thick Walls from 
Arbitrary Variation of Adiabatic Wall Temperature and Heat Transfer Coefficient, 
N.A.C.A. T.N. 4105, October 1957. 


CarsLaw, H. S. and JageGer, J. C. Heat Conduction in Solids. Oxford University Press. 
pp. 377 and 380-381, 1947. 


HEISLER, M. P. Transient Thermal Stresses in Slabs and Circular Pressure Vessels. Journal 
of Aoplied Mechanics, Vol. 20, pp. 261-269, June 1953. 


284 The Aeronautical Quarterly 


| 


Further Comments on High-Lift Testing in 
Wind Tunnels with Particular Reference to 
Jet-Blowing Models * 


S. F. J. BUTLER and J. WILLIAMS 
(Royal Aircraft Establishment, Farnborough) 


SumMMaRY: This paper discusses some of the special problems of wind-tunnel 
testing which arise with high-lift jet-blowing models, supplementing earlier 
comments’. Further detailed remarks are made about recent developments 
on tunnel-wall interference, test rigs and methods of minimising constraints 
from air-feeds to models, and on general test and model design techniques. 


1. Introduction 


Some of the special problems encountered in wind-tunnel testing of high-lift 
models were discussed by Anscombe and Williams? nearly three years ago, largely 
in terms of early experience gained by the R.A.E. and N.P.L. low-speed tunnel staffs. 
Since then, further research and development have been carried out in connection 
with tests on models incorporating plane jets (boundary layer control or jet flaps) 
and round jets. In particular, the half-model techniques already described 
previously for three-component force measurements with jet-blowing have now been 
thoroughly proved”: *:*’ and widely used. Furthermore, at R.A.E., a special tech- 
nique has been developed for six-component measurements on complete models with 
jet-blowing, and tunnel “lift-constraint” corrections have been derived theoretically 
for three-dimensional wings with full-span jet flaps. 


Further detailed comments are made on tunnel-wall interference with jet- 
blowing models, on test rigs and methods of minimising constraints from air-feeds 
to models, and on the measurement and evaluation of jet-blowing parameters. 
Earlier remarks on accuracy of manufacture requirements are supplemented by 
comments on blowing slot and duct design and on realistic representation of practical 
aircraft configurations. 


NOTATION 
a_ local velocity of sound 
A_ wing aspect ratio 


*Based on a paper prepared for the AGARD Wind-Tunnel Panel Meeting. Brussels, March 
1959. 
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A, area of slot throat or convergent nozzle exit 
b breadth of tunnel, in spanwise direction 
c wing chord 
cross-sectional area of tunnel working section 
C, jet-reaction coefficient 
see Section 4.1 
C, jet-momentum coefficient 
Cy jet-flow coefficient 
C, lift coefficient 
C,, pitching-moment coefficient 
C, thrust coefficient 
d, diameter of round-jet exit 
D diameter of circular wing 
h_ height of tunnel working section, normal to wing 
h, height of jet exit above tunnel floor 
J static jet thrust 
M, jet mass-flow rate 
Py jet (or duct) stagnation pressure 
Pp, main stream static pressure 
q, main stream dynamic head 
Re chord Reynolds number 
S wing reference area 
T, jet (or duct) stagnation temperature 
T, main stream temperature 
v; jet velocity 
V, main stream velocity 
x streamwise distance 
z distance below wing 
@ wing incidence relative to main stream 
y ratio of specific heats for air 
.o tunnel “lift-constraint” correction factors (see Ref. 9) 
¢ angle of downwash relative to main stream 
¢., angle of downwash at infinity downstream 
jet angle relative to wing chord line 
p, density of main stream 


tailplane angle 
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2. Choice of Model 


The need to minimise tunnel-wall interference, because corrections are inevit- 
ably approximate, tends to limit the size of the model relative to that of the working 
section. On the other hand, the requirements of test Reynolds number, pressure- 
plotting, and accuracy of manufacture usually demand as large a model as possible. 
Tunnel-wall interference is discussed further in Section 2.1 and a few additional 
remarks about Reynolds number effects are made in Section 2.2. Although pressure- 
plotting and accuracy of manufacture does not warrant any further comment, the 
importance of good slot and duct design and of realistic representation of practical 
aircraft configurations seems worth mention (see Section 2.3). Other conflicting 
factors which must be taken into account in model design, such as the provision of 
air-feeds to the model for blowing (and suction), the associated type and size of 
model support, the model strength, and the best use of available air supplies, are 
discussed in Section 3 of this paper and of Ref. 1. 


2.1. TUNNEL INTERFERENCE 


The various interference corrections and the corresponding limitations on model 
size were considered in Section 2.2 of Ref. 1. With boundary layer control models”, 
where the sectional lift coefficient is unlikely to be much greater than 4, interference 
considerations should be little different from those associated with normal testing. 
The ratio of wing chord to tunnel height is usually kept below one-third and the 
ratio of wing span to tunnel breadth below three-quarters. However, with models 
incorporating a large jet efflux at a high angle to the main stream, and in particular 
for tests at the higher lift coefficients associated with jet-flap schemes”, the 
interference can assume formidable proportions unless the models are much smaller 
than normal; the classical Glauert corrections’ may also be considerably in error. 
The following discussion refers to closed working sections, although some of the 
comments are of more general interest. Theoretically, partially ventilated working 
sections (slotted or perforated walls) can offer much reduced lift-constraint 
corrections and small blockage. However, the correction factors are particularly 
sensitive to the assumed effective open area, which can vary appreciably from the 
simple geometric value over a range of high lift tests. 


2.1.1. Plane-Jet Models (Jet Flaps) 


With regard to lift-constraint corrections for finite aspect ratio jet-flap models, 
the effects of induced curvature of the stream over the wing chord can probably be 
neglected as usual, but significant effects can undoubtedly arise from the alterations 
in the curvature of the jet sheet. Hence, the classical correction factors based on 
simple Glauert theory may be appreciably in error at large values of the jet 
momentum coefficient C,, and jet angle 6, when the circulation about the wing is 
controlled mainly by the effectiveness of the jet sheet rather than by wing incidence 
or camber. Maskell has recently extended the classical lift-constraint theory to 
jet-flap wings, following his treatment of the effects of finite aspect ratio for thin full- 
span jet-flap wings in an inviscid free stream™. This extended theory, for closed 
working sections, assumes that the tunnel constraint is small, and that the bulk of 


August 1960 287 


Ss. F. J. BUTLER AND J. WILLIAMS 


the bound vorticity over the wing and jet is confined within a distance from the 
tunnel centre which is small compared with the linear dimensions of the working 
section. The constraint correction can then be conveniently represented by an 
increment Az to the geometrical wing incidence 2, similar but not identical to the 
classical correction, together with an increment AC,, to the jet-momentum coefficient. 
This latter arises because of the upwash distribution induced at the centre line of the 
tunnel by the image system, causing higher curvature of the jet sheet and hence 
higher vorticity than would exist in an unrestricted stream. 


The corrections to the measured wing incidence, jet-momentum coefficient, wing 
lift and thrust coefficients, to give equivalent free stream conditions, are as follows: — 


Aa= 1 C C,, (no tail), 
AC,.= where zA+2C, ° 


AC, =AC, sin (6+ 2) +C,Az, 
AC,;=AC,, cos (06+ 2) -—C, Az. 


The measured distribution of downwash angle <(x,z) must be corrected by the 
amount 


Ac=(1+— *) Aa; 
this leads to the normal pitching-moment correction formula for tail-on results, 


x 
= 


and a correction to the measured jet path, 


Az=x(1+ Ae, 


with a corresponding correction to tail height. 


The symbols § and C denote the reference wing area and the cross-sectional 
area of the tunnel working section, respectively; 6 and 6’ represent the classical 
constraint-correction factors as tabulated in Ref. 9. It will be noted that the 
correction C,Az to lift, frequently ignored in conventional testing, can be quite 
significant with jet-flap models. In the corrections given, the use of the jet-reaction 
coefficient C, [=J/(q, S)], might be preferred to the slot momentum coefficient 
C,.[=My,2,/(q, S)], provided that the measured static thrust could reasonably be 
assumed to represent the jet momentum leaving the trailing edge with wind-on 
conditions (see Section 4.1). Furthermore, when thrust at small jet angles is being 
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@ oO 


/ 


MODEL SPAN (FT.) FOR 112 FT. TUNNEL 


10 15 


———— Aa=2°, using classical tunnel constraint corrections. 
——— Aa=2’, using new jet-flap constraint corrections, with a full-span 60° jet flap. 
Ficure |. Effect of tunnel wall interference on permissible model size. Rectangular-wing 
aspect-ratio-6 model in an 114 ft. X 84 ft. tunnel. 


studied, careful evaluation of (6+ 2-—<,) is necessary, as discussed later in this 
section. Finally, the foregoing corrections correspond strictly to the case of full- 
span jet-flap wings with an elliptic spanwise distribution of loading, so discretion 
must be used in their application to part-span jet-flap models. 


It was suggested in Ref. 1 that the model size should be limited to keep the 
incidence correction Az below 2° up to the maximum value of C, to be tested, say 
10 for S.T.O.L. research investigation. The corresponding correction, A:, to the 
angle of downwash at a norma! tailplane location might be double this value. Fig. 1 
plots the appropriate limit in size against C, (with A4z=2°) for a rectangular wing 
of aspect ratio 6 with a full-span jet flap, mounted horizontally in an R.A.E. 
114 ft. x 84 ft. closed tunnel; curves corresponding to both classical and Maskell’s 
new corrections are included. It is seen that the new corrections permit a slightly 
larger span. More important, perhaps, is the better accuracy of the new corrections, 
so that the suggested restriction on Az might be somewhat relaxed if its magnitude 
were the sole criterion. However, much broader considerations are desirable, 
particularly if an assessment of thrust efficiency forms an important part of the 
research experiments. This can conveniently be illustrated by a brief examination 
of the new lift-constraint corrections for the model and tunnel of Fig. !. 


If the span of the aspect-ratio-6 model is limited to 5 ft., to give 42 <= 2° for 
C:; <= 10, then the correction AC, to the measured C, is not likely to exceed 
0-05 C,; thus any error in this correction will not be particularly serious as regards 
lift. However, the corresponding correction —C, 42 to the measured thrust 
coefficient C; could well be as much as 0-1 C, at high C,-values, and an even 
higher proportion of C;. This could be particularly important at high C.-values and 
jet angles, when errors in the correction could give misleading answers about the net 
thrust available under S.T.O.L. conditions. Again, although the further correction 
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AC,, cos (6+ z) [=C,,42 cos (6+ z)/(6+ 2—2,.)] to the measured C; is likely to be 
less than 0-05 C,, for jet angles of 30° and above, the proportion could be higher and 
the correction more inaccurate at small jet angles, because of possible errors in the 
evaluation of the denominator (¢+2-—<,). In such cases, it is probably more 
accurate to express C,, ©, 42 and (¢+ 2—<,) all as specific linear functions of 6 
and z for the evaluation of this correction at small 6 and z. 


Thus, in limiting the model size to allow for the inevitably approximate nature 
of the lift-constraint corrections, the magnitude of Az and C, Az/C, should be 
examined for the maximum value of C, envisaged at least, together with the 
magnitude of AC,,/C,, for a representative range of values of C, and C... 


Additional interference effects can arise when the jet angle and momentum 
coefficient are so large that the jet sheet approaches close to the floor of the tunnel. 
Clearly, the foregoing lift-constraint corrections ought then to be modified because 
the jet cannot reasonably be assumed to lie along the centre line of the tunnel. The 
sensitivity of the constraint correction factors 6 and 6’ to offset from the tunnel 
centre line can be assessed using Ref. 10, but there is also a further significant effect 
associated with “jet-constraint” on the tunnel main stream flow between the lower 
surface of the wing and the tunnel floor. This usually leads to a reduction in the 
measured lift for a particular configuration, in contrast to normal lift-constraint 
effects; pitching moment, downwash, and drag changes also become embarrassing. 
As suggested in Ref. 1, effects should be small with a wing chord not exceeding 
one-sixth the tunnel height, at least for jet angles 6 up to 60° with C,,-values up to 
4 (Ref. 11). At larger 6 and C.,, say up to 90° and 10 respectively, a wing chord as 
little as one-tenth the tunnel height could be needed“**. Fortunately, for our 
deliberations, the lift-constraint rather than jet-constraint limitations are usually the 
more significant in practical cases. For example, the representative aspect-ratio-6 
model in the 114 ft. x 84 ft. tunnel is limited by lift-constraint considerations to 60 in. 
span, so that the 10 in. chord is only one-tenth the tunnel height. More guidance 
on the effect of proximity of the jet sheet to the tunnel floor should be obtained from 
current R.A.E. tests on a jet-flap complete model with various ground heights. 


2.1.2. Round-Jet Models 


Tunnel interference corrections for models with round jets at large angles to 
the main stream have so far received little attention, since the problem is a fairly 
novel and difficult one. If the jets are of sufficient spanwise extent to produce 
substantial jet-flap effects, then the foregoing arguments could of course be applied 
qualitatively. In many practical cases, however, the pressure lift generated on the ex- 
ternal wing surface due to the jet efflux is only a fraction of, rather than many times, 
the gross jet thrust and may even be negative. The path of the jet is then primarily 
a function of the ratio 7,/V, of jet to tunnel speed and is not intimately associated 
with increased bound vorticity on the wing, as in the jet-flap case. The tunnel 
interference is therefore mainly associated with the round jets hitting or approaching 


*Mounting the model above the tunnel centre line should alleviate the jet-constraint effect 
but would introduce further complications and unknowns. 
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close to the floor, so that conventional lift-constraint arguments become quite 
inappropriate. In the absence of theoretical guidance, test results on ground effect 
could again be interpreted to suggest limitations on model size, but unfortunately 
the available information is almost all restricted to zero forward speed. 


With wing areas only slightly larger than the jet exit area, the ground inter- 
ference at zero forward speed can be mainly attributed to back pressure on the jet 
and the lower wing surface. For most test purposes, this effect is probably negligible 
when the jet exit diameter d [= /(4A,./=)] is less than one-half the height A, of the 
jet exit above the floor, i.e. less than one-quarter the tunnel height for a centrally- 
mounted model. With gross wing areas several times the jet exit area, the 
interference at zero forward speed largely arises from suctions on the wing lower 
surface associated with flows induced by entrainment with the jet. In R.A.E. tests 
on wings of circular and equilateral plan forms, having a centrally mounted jet with 
a ratio D/d of wing to jet diameter exceeding 4, Wyatt has shown that the ratio of 
the downward suction S, on the wing surface to the gross thrust 7 of the jet is 
mainly a function of D/h;. The ratio S,/T decreased from about 0-2 to 0-06 as D/ h, 
was reduced from 4 to 1:6, and probably would have been less than 0-02 with 
D/h,<1/3. With annular rather than circular jets, the interference problem 
becomes somewhat more complex, since downward suction on the peripheral wing 
surface and back pressure on the central fairing or hub can be present simul- 
taneously. Again, if air intakes are represented on the jet-blowing model, they may 
be subject to interference due to re-circulation of the jet efflux and, in this respect at 
least, the roof of the tunnel can also play an important part. 


When the main stream speed is non-zero, tunnel interference effects on net thrust 
and pitching moment will be equally as important as those on lift. In the absence 
of adequate information, it seems reasonable to conjecture that serious interference 
is unlikely for most practical cases if the wing chord does not exceed one-third of the 
height of the jet exit above the tunnel floor, i.e. one-sixth of the tunnel height. This 
limitation satisfies the jet-constraint restriction associated with jet-flap arguments, 
should meet most lift-constraint restrictions associated with pressure lift (and drag) 
on the wing arising from partial jet-flap effects or wing incidence, and should give 
stai lifts correct to within one or two per cent. Admittedly, it may be rather severe 
for many model tests with round jets and often may have to be relaxed, particularly 
for wing-fan models. Each case has naturally to be considered on its own merits, 
both in respect of the information needed and the accuracy required. Check tests 
can be made by varying the height of a ground-board but, as stressed in Section 3.3 
of Ref. 1, care is needed to ensure that the different blockages in the two compart- 
ments of the tunnel do not lead to misleading results. 


To provide further information on tunnel interference with round-jet models, 
some tunnel tests are being made at the R.A.E. on a simple single-jet model at 
various heights above the tunnel floor and with various ground-board heights. In 
addition, some of the more complex powered-lift models now being used for basic 
research in the 114 ft. x 84 ft. and 13 ft.x9 ft. tunnels at R.A.E. are also being 
briefly tested in the 24 ft. diameter open-jet tunnel to give a check on tunnel 
interference effects. 
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2.1.3. “Blockage” Corrections 


In high-lift testing, solid and wake blockage corrections to the tunnel dynamic 
head are usually only of secondary importance when choosing model scale. Never- 
theless, they should not be overlooked, particularly if good absolute accuracy is 
required in, say, the determination of C,.m.x or the thrust coefficient. Solid blockage, 
associated with the volume of model in the bounded main stream, can be readily 
estimated using classical methods, but “wake blockage” can arise in severa! novel 
forms with jet-blowing high-lift models. 


When the flow is unseparated on conventional models, the wake blockage 
correction to the main stream dynamic head q, can conveniently be expressed in 
terms of the rig drag coefficient Cp, and the model profile drag coefficient C»,, 
as follows :— 

Ago! =4 (S/C) (Cor + Cov). 


For completely stalled flows, unpublished work by Maskell at the R.A.E. 
suggests the relation 


Aq, / 4 (S/C) Co. )+ 1K (S/C) (Cy 


where C,,, again represents the profile drag coefficient with unseparated flow, Cy, 
represents the induced drag coefficient (C,*/=A), Cy denotes the measured total drag 
coefficient under the relevant test conditions, and K is an empirical constant of about 
5 for aspect ratios up to 10. For intermediate cases, where some degree of 
separation is present although the flow is not completely stalled, the choice of a 
K-value (between 5 and 1) is not a straightforward one and will depend on whether 
the separation is of the trailing-edge or nose-bubble type. 


On jet-blowing models, apart from such low velocity wakes, consideration is 
also necessary of the high-velocity jets which are being fed into the tunnei main 
stream. These jets steadily grow in width and decrease in maximum velocity as 
mixing takes place with the lower velocity main stream. In the particular case of 
jet flaps, any separation on the wing upper surface usually reattaches forward of the 
trailing edge and the thickened boundary layer is completely absorbed by the jet 
emerging close to or at the trailing edge; the resulting “thick wake” is quite 
different in character from that behind a conventional aerofoil. The problem of 
blockage corrections with jet-blowing models would therefore appear to be rather 
complex and to warrant further investigation. 


2.2. REYNOLDS NUMBER 


The earlier remarks“? concerning the requirements of test Reynolds number 
have been well substantiated. In the absence of appreciable flow separation on the 
wing upper surface, variation of wing Reynolds number has had no significant effect 
on jet-blowing wing performance over the range tested (R.=0°5 x 10° to 2x 10°). 
The main stream Mach number, although low, seems likely to be more important, 
since sonic flow can occur over the wing surface at high lifts. The tailplane (or fin), 
on the other hand, may be working at extremely low Reynolds numbers and over a 


292 The Aeronautical Quarterly 


HIGH-LIFT WIND-TUNNEL TESTING 


wide range of downwash angles, if fixed, so it is probably subject to considerable 
scale effect and should be used only as a mean downwash (or sidewash) meter. The 
full-scale tailplane (or fin) characteristics should then be estimated using full-scale 
lift slopes and stalling angles. 


It is worth adding that the maximum possible test Reynolds number for any 
prescribed value of the jet-momentum coefficient C,, is independent of the model 
scale, if the available air supply is being used to its limit. Likewise, the maximum 
possible lift and thrust forces on the model are also independent of scale, if C,, is 
the parameter governing the effectiveness of blowing on a specific configuration; 
however, the corresponding moments will be proportional to model scale. 
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FiGure 2. Typical jet-blowing slot arrangements. 


2.3. MODEL DESIGN 


Simple convergent blowing slots have worked well over a wide range of pressure 
ratios and, with most models, the incorporation of complex convergent-divergent 
slots for supercritical pressure ratios is probably unnecessary. Some typical slot 
designs are shown in Fig. 2; Tufnol and Araldite have often proved more convenient 
materials than steel. Usually it is convenient to have a short parallel length of slot 
at the exit so as to ensure that the jet is directed properly, and also to accommodate 
the simple spacers needed to maintain the correct slot width. Much larger parallel 
lengths than about five slot widths are undesirable in view of the losses due to the 
slot boundary layer. The spacers should be as small and as widely-spaced as 
practicable, consistent with maintaining the desired slot width; say, not more than 
0-25 in. wide at spanwise intervals not less than 3 in., giving an overall blockage less 
than 10 per cent, particularly for boundary layer control at large flap angles. The 
slot upper lip can advantageously be kept thin and ridges on the flap upper surface 
should be kept as small as possible. Simple cylindrical ducts are adequate to 
maintain a satisfactory spanwise distribution of slot velocity if the chordwise cross 
section of the duct is large enough to keep the duct spanwise velocity below about 
two-tenths of the slot velocity. Where this is not possible, careful tapering of the 
duct is required to ensure a uniform spanwise distribution of duct static pressure, 
and corner vanes are needed to turn the air chordwise if total loss of duct dynamic 
head on turning cannot be tolerated. 
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In testing models of specific aircraft configurations, the effect of essential 
practical features should be examined; for example, flap gaps, flap cut-outs, and 
extensive blockage of the slot such as might arise from the flap operating mechanism. 
The effects of air brakes, undercarriage assemblies and stores can also vary 
significantly with jet-blowing. 


3. Various Test Rigs and Methods of Air-Feed 
3.1. HALF-MODELS 


The advantages of using a half-model are that a larger-scale model can be 
installed, that the manufacturing time is considerably reduced, and that the problems 
of model support, method of air-feed, and carriage of pressure-plotting tubes are all 
greatly simplified. The disadvantages are that the model cannot be yawed, so that 
the measurements are therefore restricted to lift, drag, pitching moment, and possibly 
rolling moment, at zero angle of yaw. 


The two half-model techniques described in Section 3.1 of Ref. | have now been 
used successfully for a variety of blowing tests, including jet-flap models. Basically, 
the N.P.L. method (Fig. 3) is designed for a simple three-component balance, with 
the model mounted horizontally on struts or wires, and a vertical false wall inserted 
in the tunnel. The horizontal air-feed is aligned with the model pitch axis, being 
designed to minimise air-feed constraint effects on lift, drag, and pitching moment. 
Variants of this method have been used for similar tests in several other wind tunnels 
in Great Britain, with equal success. The R.A.E. method (Fig. 4) is suitable for a 
virtual-centre balance, preferably a floor balance, with the model mounted rigidly 
from the moment table*, using the tunnel floor as a reflection plane. The vertical 
air-feed on the model pitch axis is again designed to minimise air-feed constraints 
on the relevant force measurements. 


*With a virtual-centre balance, the model can of course be connected to the moment table 
in any convenient way, frictionless mountings being unnecessary. This is especially advantage- 
ous with air-feeds to the models. 
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As regards air-feed connectors for half models, the “necked canvas-hose” and 
“two-faced air-bearing” types developed at R.A.E. and N.P.L. respectively have 
already been described in Ref. 1. Some other types are mentioned in the general 
discussion on air-feed connectors in Appendix I. 


The results obtained from half-model tests have shown good agreement with the 
corresponding results obtained on complete models. However, comparative high- 
lift experiments on a half-model and a similar complete model in the same tunnel 
are still lacking. From the broad comparisons available, which include in one case 
a jet-blowing configuration, it would appear that the tail-off lift, drag, and pitching 
moments are generally in good agreement, although significant differences can arise 
in zero-lift drag and pitching moment. The downwash « at the tail and d:/dz, as 
determined by using the tailplane as a mean downwash meter, agree within experi- 
mental accuracy. However, as might be expected, tail power may be different, 
because of wall boundary-layer effects on the half-model, and possibly to scale effect. 
In any case, as discussed in Section 2, it is advisable to use the tail mainly as a 
downwash meter. 


The reflection-plate boundary layer is normally of little concern, provided that 
the half-model has a fuselage thickness rather more than three times that of the 
boundary layer, and provided that the static pressure gradient in the tunnel is not 
sensibly affected. For tests on a half-wing without a fuselage, it may be necessary to 
thin the reflection-plate boundary layer, at least near the wing root. This could be 
achieved by using suction to remove the boundary layer or a vortex-generator 
system to deflect it. 


Clearly, the half-model technique can yield valuable results in less time than the 
complete model. However, to enable lateral stability in sideslip and asymmetrical 
aircraft configurations to be studied, it is also essential to dev-lop complete-model 
rigs suitable for six-component measurements on jet-blowing models. 


3.2. COMPLETE MODELS 


The additional difficulties which arise in designing a complete-model rig with 
jet blowing are quite appreciable, particularly as regards achieving completely 
satisfactory model support and air-feed arrangements. Some possible solutions 
based on recent R.A.E. experience are therefore discussed. 


3.2.1. R.A.E. Jet-Blowing Complete-Model Rig and Test Technique 


The jet-flap complete-model rig briefly mentioned at the end of Section 3.2 of 
Ref. 1* has now been fully developed and in satisfactory operation for more than 
two years; the rig also permits tests on floor-mounted half-models over a 40° 
incidence range. Complete models are mounted horizontally (Fig. 5) on a 
single hollow strut from the virtual-centre lower balance of the No. 2 114 ft. x 
84 ft. tunnel, the air for jet-blowing being fed to the model through a “four-faced 


*Mr. A. Anscombe was responsible for the initial design of the rig. 
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Ficure 5. General arrangement of R.A.E. complete-model rig. 


air-bearing” connector specially devised to minimise the constraint effects on all six 
components. Details of this air-bearing connector are given in Appendix II. The 
tunnel balance, briefly described in Chapter V of Ref. 9, is also unusual in that it, 
too, is supported on compressed-air bearings. 


The main features of the complete-model rig, apart from the air-bearing 
connector, are concerned with the pitching and yawing mechanisms at the centre of 
the model, frictionless bearings being unnecessary since a virtual-centre balance is 
being usedt. 


The model is pivoted on a hollow central distributor (Fig. 5) and can be pitched 
relative to the distributor. To ensure rotatable leak-proof connections, O-ring seals 
are incorporated between the twin output pipes of the distributor and the main wing 
ducts which rotate with the model. The distributor is bolted by its base flange to 
the top flange of the hollow support strut, slotted holes in the flange permitting up to 
+ 15° yaw of model and distributor relative to the strut; a dowel can be used to set 
relative yaw at 0°. An O-ring seal is again used to avoid leaks. The strut access 
hole in the bottom of the fuselage is sealed by a roller blind to avoid possible 
“ bomb bay ” flow interference effects. The top of the strut, which is exposed during 
longitudinal stability tests, is faired to a low-drag 40 per cent aerofoil section and the 
strut is mounted on a main base plate spanning the balance pedestals, which have 
been placed side by side for convenience. The strut guard is mounted independently 
on the floor turntable. 


+See footnote on p. 294. 
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The air-bearing connector is isolated from the model and strut by the main 
base plate, so that strut deflections under load are not transmitted directly to the 
balance side of the connector; this occupies the space below the base plate and 
between the pedestals. Appendix II contains details of the final arrangement (see 
Figs. 6 and 7). 


During longitudinal stability tests over an incidence range at zero sideslip, the 
model incidence is controlled by a tail wire connected to a capstan mounted on the 
balance, with a loaded front wire to return the model towards negative incidence. 
The results are corrected for drag and pitching-moment tares due to the exposed 
portion of the strut and the wires. 


During lateral stability tests, the model incidence is fixed by a dowel-pin and 
the front wire and tail wire are removed. It has been found preferable, with current 
jet-flap models, to extend the strut guard so as to cover the whole of the exposed 
strut, and to rotate the model and strut together by yawing the moment table*. In 


*At the model design stage, it was not known whether the moment table could be yawed 
without introducing serious air-feed constraint effects on one or more of the freedoms; so the 
internal yawing joint was included. 
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this way, strut tare effects on side force and rolling moment are avoided, and the 
moments are measured about yawed axes directly. Moreover, the balance zeros do 
not vary greatly with yaw angle, as would be the case if the model were yawed 
relative to the strut. Although the interference of the extended strut guard on the 
fin and tail must be larger than that of the exposed strut, it has been shown to 
be small. 


The position of the strut guard giving least interference, as indicated by pressure 
plotting on the strut, is virtually along wind. Unfortunately, the size of the strut 
guard permits only 7° yaw relative to the strut, so that at large model yaw angles 
the strut guard is somewhat misaligned from the optimum. By comparative tests at 
a fixed model yaw angle, with a range of strut guard yaw angles, the measured 
results can be extrapolated to the zero strut guard yaw condition. For the present 
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Ficure 8. Possible alternative to air-bearing system of Fig. 7, using simple flexible connectors. 


model, the errors due to misalignment were relatively insignificant, even with the 
model at incidence and with the fin and tail on. 


The tunnel results are corrected in all cases for the calibrated effects of shift 
of the model centre of rotation relative to the balance virtual-centre, since the model 
support strut can deflect under load. On the present model the deflections involved 
are small, generally less than 0-020 in. Another small correction is applied for the 
calibrated misalignment of the model centre of rotation and balance virtual-centre in 


rigging. Balance zeros are taken with the correct static pressure in the air-feed 
connector, the supply being sealed off from the model blowing ducts by an 
isolating valve. 


Since the zeros do not change rapidly with internal static pressure for the final 
system, the resulting zeros are only slightly less accurate than those for the basic 
balance. It has also been established that the balance calibration is not affected 
appreciably by the presence of the connector, and that the constraint effects due to 
the connector remain the same when loads are applied to the model. Such questions 
would, of course, require examination with any connector system; it cannot be tacitly 
assumed that the presence of the air connector will merely affect the basic zeros. 


3.2.2. Possible Extensions of the R.A.E. Complete-Model Rig to other Balance 
Systems 


Since the R.A.E. complete-model rig has proved very satisfactory in practice, 
some possible modifications to permit six-component measurements for jet-blowing 
models on other balance systems seemed worth mention. 


With a virtual-centre balance, either above or below the tunnel, the R.A.E. rig 
should require little modification, provided the moment table can be assumed to 
return to sensibly the same position and attitude when balanced under typical aero- 
one further degree of rotational freedom in pitch, as discussed later, in order to 
dynamic loads as when balanced at zero load. It may be necessary otherwise to add 
permit slight misalignments of the “balance side” relative to “earth side” in each of 
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Ficure 9. Some possible air-feed arrangements for use with conventional (not virtual-centre) 
six-component balances. 


the three translational and three rotational senses without introducing substantial 
constraints. Again, when the table displacements are large, inside and outside 
diameters of opposing faces of the air-bearing connector may have to be made 
unequal. The air-bearing connector could probably be replaced by an arrangement 
using three flexible connectors, such as bellows or necked canvas sleeves, as shown 
in Fig. 8, provided larger constraints were acceptable; a turntable would have to be 
incorporated on the “balance side” or “earth side” to permit the balance to be 
yawed. This system could also completely avoid the air leak present with 
air-bearing connectors or labyrinth seals*. 


With a conventional balance (not the virtual-centre type), the air must be fed 
directly to the model through a frictionless connector. Furthermore, greater move- 
ments of the model under load will usually occur, both in translation and in rotation. 
Some possible arangements using air-bearing connectors are outlined in Fig. 9, for 
models supported from an overhead balance in the conventional manner. 


The first scheme employs a hollow strut, suspended from the model by a 
spherical air-bearing at the model centre of rotation. A “four-faced air-bearing” 
connector (or bellows system) is used at the foot of the strut, with unequal inner and 
outer diameters on opposing faces to permit substantial relative translational 
movements without appreciable alteration in force constraints. Thus, as the model 
moves, the strut will remain vertical, with its weight acting through the model centre 
of rotation, while the spherical bearing allows for model rotational movements, 
including those due to model distortion under aerodynamic loads. The main 


*With air-bearing connectors, the leak could be prevented, if desired for mass flow calibration, 
by incorporating inflatable rubber seals in the plane faces of the air-bearings. 
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objection is the aerodynamic interference which could arise due to the air leak from 
the spherical bearing inside the model; the mass flow calibration in itself should 
cause no difficulty since it could be carried out on a special rig or after sealing the 
air gaps. The second scheme employs the same lower part, but the spherical 
bearing inside the model is replaced by a cylindrical joint and O-seals, as in the 
R.A.E. complete model. Preferably, the top of the strut should be maintained 
vertical automatically as the model incidence is altered. Bellows or a necked-canvas 
joint would be incorporated as near the top as possible, so that the lower part of the 
strut would be able to align itself truly vertically with its c.g. practically below the 
model centre of rotation, even if the top of the strut were slightly misaligned in pitch 
or roll from the vertical. The third scheme dispenses with the flexible joint, but 
might require the additional rotational freedom in pitch, if the lower part of the strut 
cannot be aligned sufficiently accurately: 


All these schemes appear rather complicated, but would at least be expected to 
result in relatively small air-feed constraints on the model, and for that reason would 
prove preferable to cruder schemes involving larger constraints. The advantages of 
virtual-centre balances for six-component measurements on jet-blowing models are 
obvious in view of the increased complexity otherwise involved in the air-feed 
connector system. 


3.3. Two-DIMENSIONAL AEROFOIL MODELS 


Investigation of two-dimensional characteristics, to check basic theory and to 
study specific fundamental aspects, usually requires overall force measurements as 
well as pressure plotting and wake traverses. The characteristics could be derived 
from results obtained on a finite aspect ratio model, but the analysis is likely to be 
difficult and unreliable unless the aerofoil is mounted with reflection plates at its 
ends. Some of the problems which arise in testing such a model between reflection 
plates at high lift coefficients, have already been discussed in Section 3 of Ref. 1. If 
the model completely spans the tunnel from wall to wall, then wall boundary-layer 
separations and flow through the air gaps at the ends of the aerofoil are two possible 
sources of trouble. When it is impracticable to span the tunnel, because of 
limitations of model strength or air supply, separate end-plates can be used as false 
walls bridging the tunnel from roof to floor (assuming a horizontal model). 
Difficulties then arise in ensuring a uniform flow through the central compartment 
acting as a “new working section,” and in determining the correct main stream 
velocity*. If the end-plates are much smaller, so that comparting the tunnel is 
avoided, the aerofoil rig becomes only quasi-two-dimensional in the sense that the 
chordwise loading is sensibly independent of spanwise position but downwash effects 
are still present (see Appendix IV of Ref. 1). 


A better test arrangement comprises an aerofoil spanning the smaller dimension 
of the tunnel, with the centre section of the span mounted and fed independently of 


*These difficulties could be avoided by fitting a complete “ two-dimensional ” contraction and 
diffuser to the false walls so that the whole of the tunnel main stream flow passes through 
the new working section. But special care is then needed to avoid flow separation in the 
diffuser and deterioration in the tunnel fan performance. 
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FiGure 10. Possible two-dimensional test arrangements for R.A.E. virtual-centre floor balance. 


the twin outboard (dummy) sections, but with the spanwise distribution of aerofoil 
shape and jet blowing as uniform as possible. The small gaps between the centre 
and outboard sections should be sealed, without introducing significant constraints 
when balance measurements are being taken (Section 3.3.2 of Ref. 1). Moreover, 
the associated breaks in the spanwise extent of the blowing slot should be kept as 
small as possible, about the same size as the width of spacer being used to set the 
slot gap. Since the outboard sections are intended to avoid tunnel-wall boundary- 
layer effects on the centre section, the geometrical aspect ratio of the whole aerofoil 
must be kept as large as possible (say, at least 2), commensurate with the need for a 
large chord, because of Reynolds number, pressure plotting, and manufacturing 
considerations. In any case, flow separation at the tunnel walls should be minimised 
by auxiliary boundary-layer control on the walls themselves. The chief difficulties 
likely to arise with this type of two-dimensional rig are the support of the centre 
section independent of the outboard ones, the separate air-feed to the centre section, 
and the accurate alignment of the centre and outboard sections. These difficulties 
can be somewhat eased by keeping the geometrical aspect ratio of the whole aerofoil 
low, commensurate with the needs already mentioned and with limitations of 
chord / height ratio because of tunnel interference. 


Some of the rigs designed for two-dimensional jet flap tests in the R.A.E. No. 2 
114 ft. x 84 ft. tunnel are now briefly outlined, by way of illustration. In a rig 
designed and built about four years ago, the model span had to be limited because of 
air-feed and strength considerations, and end-plates spanning the tunnel were used. 
The arrangement proved unsatisfactory for accurate measurements, largely because 
the size and blockage of the outer tunnel compartments were considerably different 
from those of the central one containing the jet-flap aerofoil. Since then no further 
balance measurements have been made on two-dimensional models. But possible 
new rigs spanning the tunnel have been considered, taking advantage of the higher 
air-feed pressure ratios which should soon be available (up to 75:1). The 
geometrical aspect ratio of the aerofoil needs to be at least 6b/h, to meet the 
limitation c/h < 1/6 of Section 2.1.1. 


The horizontal-model rig shown in Fig. 10(a) closely approximates to the ideal 
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arrangement just discussed, except in that A is the smaller dimension (84 ft.) of the 
tunnel. The centre test wing could be supported by two hollow metal struts, suitably 
covered by strut guards; the wing would be able to pitch, using a similar arrange- 
ment of metal bearings and O-seals as for the R.A.E. complete model, with the air 
supplied through the struts. Large total head losses would be involved, and further 
losses would be introduced if throttling of the model blowing duct were essential to 
obtain a uniform spanwise distribution. The spanwise location of the struts is 
chosen to maintain the extremities of the test rig horizontal with uniform spanwise 
loading. The calculated deflections under load are very small (less than 0-01 in. at 
C.=5, V,,=100 ft./sec.), and it should not be necessary to adjust the vertical height 
of the outer wings. The latter would be mounted on auxiliary turntables with 
independent air supplies and inflatable rubber seals could be incorporated at the 
joints of the wing and flap contours. 


The vertical-model rig shown in Fig. 10(b) is a simple extension of the normal 
half-model arrangement, and makes h the larger tunnel dimension (114 ft.). How- 
ever, the deflections would be embarrassingly large, about 2 in. at the tip at 
C.=5, V, =100 ft./sec., unless a restraining load were applied. Since both ends of 
the model are located in the tunnel boundary layer, repeat balance measurements 
with the model reduced in span and one dummy wing would be desirable, to 
determine the appropriate corrections. 


4. Measurement of Jet-Flow Parameters 


The main parameters and their experimental determination were discussed fully 
in Section 4 of Ref. 1, formulae for the jet-momentum coefficient and notes on 
accuracy of measurement being also included in Appendix V of Ref. 1. Some 
qualifications of the earlier comments are made here in the light of more recent 
experience. 


4.1. DEFINITION OF THE MAIN PARAMETERS 


The jet-momentum coefficient C, [=M,,/(q,S)], rather than the jet-quantity 
coefficient C, [=M,/(p,V.,S)] or the velocity ratio v,/V,, has been established as a 
satisfactory non-dimensional parameter for the correlation of basic studies on plane 
jets with arbitrary slot width*. This procedure has been justified for both subsonic 
and supersonic blowing, even up to pressure ratios as high as 5 and with simple 
convergent nozzles. However, some elaboration may be necessary when the pressure 
ratio is so low that the jet velocity is only about twice the main stream velocity V, 
and barely exceeds the local stream velocity V near the slot exit. For example, in 
boundary layer control studies, the rate of mixing is important and the relative 
velocity 7, —V becomes more significant than v, itself; in fact low-pressure boundary 
layer control results tend to correlate better with those at high pressure ratio if 
we take 

C.=M, (v7; V)/(q, S)=C. V/V). 


*For basic studies, mean sectional coefficients are usually derived by taking the wing reference 

area as that corresponding to the spanwise extent of the blowing slot. For performance 
studies, comparing differing extents of blowing, overall coefficients may be preferable, based 
on the same gross wing area as the overall force coefficients. 
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Ficure 11. Theoretical variation of jet-blowing parameters with pressure ratio. 


Unfortunately, the difficulty arises in specifying the values of V/v, or V/V,, 
particularly since V can vary rapidly chordwise near the slot exit or over the 
flap knee. 


For the evaluation of C,, the blowing velocity 7, is perforce determined 
theoretically, by assuming isentropic expansion from the total pressure p, (and 
temperature 7,,) in the slot throat or duct to the main stream static pressure p,. 
Total head variation across the width of the slot due to slot boundary layer growth 
is not usually significant from high-lift aspects if the slot is well designed (see 
Section 2.3), but should not be overlooked when a precise assessment of thrust losses 
has to be undertaken. The mass flow rates can likewise be estimated theoretically 
for the prescribed pressures, slot temperature, and slot throat area A,. However, in 
experimental work, the mass flow is best taken as that actually fed to the model duct, 
so that errors due to throat area determination and viscous effects are avoided. 


For completeness, Fig. 11 illustrates the variation of 7,/V, with blowing 
pressure ratio py»/p,, at representative conditions of main stream speed 
ft./sec.), ratio of slot temperature to main stream temperature (7/7, = 1) 
and the speed of sound in the main stream (a, = 1,117 ft./sec.), taking y=1-4. The 
corresponding values of C, and Cy are also shown for a representative ratio of slot 
area to wing area (A,/S=5 x 10~*) and the procedure for scaling to other conditions 
is indicated. 


The jet-reaction coefficient C, [=J/(q,S)] is sometimes used instead of C,, 
where J represents the total jet reaction measured wind-off with a correction added 
for the measured pressure drag round the model arising from the local flows induced 
by the jet stream. With trailing-edge slots, or simple round nozzles, C, should be 
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only two or three per cent different from C,,, but the difference may be appreciable 
if the air is ejected over a trailing-edge flap, because of skin friction losses. The 
value of C, yields the installed static thrust and, subject to some qualifications, 
signifies the effective jet momentum leaving the trailing edge. For example, the flap- 
drag correction can be quite appreciable (about 0-05 C,,) and its determination on a 
three-dimensional model is not easy, particularly if a body is present. Furthermore, 
the momentum leaving the trailing edge wind-on is not necessarily the same as that 
wind-off, for a prescribed practical condition. The value of C,, yields, by definition, 
the mass flow and total pressure required at the slot with wind-on conditions and, 
subject to the specification of w,, signifies the momentum ejected from the slot. 


It is commonly argued that, provided the jet-momentum or reaction coefficient 
is correctly represented, temperature differences between the jet and main stream 
are of secondary importance as regards model force measurements. While this may 
well be true over the temperature range of 15°C to 50°C normally practicable in 
tunnel testing without special precautions, the extrapolation of such test results to 
full-scale jet temperatures certainly needs further justification. In fact, the whole 
question of the interpretation of test results at very low forward speeds warrants 
careful examination when the jet velocity and temperature are substantially different 
from full-scale values. 


4.2. EXPERIMENTAL DETERMINATION OF MOMENTUM COEFFICIENT 


Measurements have to be made of the total head in the slot or nozzle (for the 
evaluation of v,) and of the rate of mass flow M, to the model, as discussed in 
Section 4.2 of Ref. 1. The spanwise variation of slot total head can be determined 
accurately by traversing a flattened Pitot located in the parallel slot exit, provided 
the slot is wide enough for the purpose, say, 0-01 in. or more. At the same time, the 
spanwise variation of slot width can usefully be checked with feeler gauges, prefer- 
ably over the same range of jet pressure and temperature as in the tunnel tests. A 
static pressure tapping in the low velocity duct can be used as a convenient datum 
for setting the nominal pressure ratio. If the slot width is too small to admit entry 
of a flattened Pitot, the mean pressure ratio, as indicated by several such pressure 
tappings, has to be used. In this case, it is especially important that the duct 
velocity should be low enough to ensure negligible losses and a uniform spanwise 
distribution of slot total head. 


Although the mass flow can be readily measured by an orifice plate in the 
delivery line between the pump and the model, its determination to within an 
absolute accuracy of +1 per cent over a wide range of test conditions is no easy 
matter. Fortunately, a lower absolute accuracy suffices for most purposes, provided 
the comparative accuracy can be maintained at this figure (see Appendix V of 
Ref. 1). The precise assessment of thrust losses, particularly under cruising 
conditions, represents one case where a good absolute accuracy in the measurement 
of M, and C,, may be essential. 


If the mass flow meter has to be installed in non-standard conditions, as often 
happens because of space limitations, then it must be recalibrated in situ and any 
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leaks downstream of the meter must be sealed or allowed for. The measured thrust 
characteristics of a simple circular nozzle of about the same exit area also form a 
useful standard of reference against which to compare My,. Finally, it is advisable 
to measure the mass flows for the prescribed duct reference pressures under 
completely representative model and wind-on conditions. For example, with sub- 
sonic blowing over flaps at large deflection angles, the mass flow can be noticeably 
affected by the variation of the static pressure just outside the slot with the main 
stream velocity and flap angle. 


Appendix I 
Some GENERAL REMARKS ON THE CHOICE OF AlR-FEED CONNECTORS 


In view of the variety of possible air-feed connectors, the merits of some known 
to the authors are discussed briefly here. Clearly, no single type will provide a universal 
solution for all model tests, but the following remarks may assist in the selection of 
a suitable type for a particular application. 


A mercury seal is attractive when sufficient vertical height is available, since this 
limits the blowing pressure ratio which can be tolerated. Difficulties may arise from 
the weight of mercury involved and the problem of maintaining the strut in a truly 
vertical position. Significant tare forces could also occur from the internal static 
pressure, and from vertical momentum and static pressure changes along the strut. associ- 
ated with surface friction and separation, if the velocity is high. A simple rubber tube 
hanging in a U-shape may likewise be adequate if sufficient vertical depth is available, 
particularly when the air can be supplied at high pressure so that the tube diameter can 
be small. Measurable constraints seem inevitable in this case. 


A quite different approach is to use a hollow strain-gauge sting for the model 
support. Here the principal objections are the necessity for suitable electronic equipment, 
the problem of adequate accuracy over a wide range, including difficulties associated 
with heating effects of the air supply, and possibly the lengthy development time. 


A multiple arrangement of bellows and gimbal rings would seem attractive, but 
the authors have no personal experience as to the magnitude of the resulting constraints, 
or as to whether substantial model movements could be tolerated. The simple “* necked- 
canvas sleeves,” mentioned in Appendix II of Ref. 1, are essentially flexible bellows 
and could be employed along the lines suggested in Section 3.2.2 to permit six-compon- 
ent force measurements. Appreciable constraints could arise, but with good design they 
should be repeatable and could be allowed for. The resulting system should be robust 
as well as simple and probably quite adequate for most purposes. 


Air-bearing connectors are discussed in detail in Appendix III of Ref. 1 and 
Appendix II of this paper. A possible disadvantage, compared with bellows or hose-type 
connectors, arises from the inherent air leak, but this need not be troublesome, except 
perhaps when the air-bearing has to be located inside the model itself. Although appear- 
ing somewhat complex at first sight, the parts required are reasonably straightforward to 
make with a lathe and shaper, provided that flat faces only are needed. Those developed 
to date at the R.A.E. and N.P.L. have proved quite satisfactory and well worth the extra 
effort required. The likely advantages comprise trouble-fee operation if the air supply is 
filtered, small and repeatable air-feed constraints to the model, and compactness. 
Furthermore, the air-bearing connector can sometimes be designed to incorporate the 
required rotations of the balance side of the connector relative to the earth side arising 
when the model is yawed or pitched, thus avoiding the necessity for an auxiliary 
turntable 
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Appendix II 
DETAILS OF “ Four-FacepD A1R-BEARING CONNECTOR FOR R.A.E. COMPLETE-MODEL RIG 
(Suitable for six-component measurements using a null-movement virtual-centre balance) 


An air-bearing connector suitable for three-component tests had already been 
developed at N.P.L. (see Appendix III of Ref. 1). In the R.A.E. jet-flap complete-model 
tig (see Fig. 5) the air-feed connection had to be chosen to minimise the constraint 
effects on all six components. Small constraints were acceptable, as long as they were 
repeatable and measurable. Since it was desired to continue testing half-models on 
the same balance, the chosen system had to permit the balance to yaw over a total 
range of about 40°, without serious constraint variations for half-model tests. Two 
main schemes were tried in turn (see Figs. 6 and 7), each using air-bearing connectors. 


In Scheme A (Fig. 6) two horizontal air-bearings were used, as in the N.P.L. three- 
component rig. The “earth side” supported a hollow Dural central block, which was 
arranged with its top and bottom faces truly horizontal. The “ balance side” faces 
of the air-bearing were formed by two machined circular pads, spaced apart by three 
shouldered round bars and hung from the base plate by three like bars. 


However, it was found that, with this scheme, the lift-zero varied rapidly with 
differential gap (G,—G,) at a fixed internal pressure. For example, a change of 0-001 in. 
in (G,—G,) altered the lift constraint by 10 lb. at an internal pressure of two atmos- 
pheres (gauge). The other constraints did not vary with (G,—G,) and were solely 
functions of internal pressure; they were large but not intolerable, as seen from the 
following table. 


BALANCE ZERO CONSTRAINTS AT INTERNAL PRESSURE OF TWO ATMOSPHERES GAUGE 


Side Pitchine Rolling Yawine 
Lift Drag force moment moment moment 
(Ib.) (Ib.) (Ib.) (Ib. ft.) (Ib. ft.) (Ib. ft.) 


Veriable <0°5 <0°5 <10 <§ 
<0°5 <0°5 <0°5 <2 <4 <0 


In view of the indeterminate lift-zero constraint with Scheme A, Scheme B was 
next tried. Here, two vertical air-bearings were added (Fig. 6) so that the block now 
floated freely between the balance side and the earth side. As in the first scheme, 
the air-bearings were supplied by direct leak from the main internal feed; however, for 
convenience, two rings of holes in the bottom pad were supplied from the interior to 
start the block floating. From the table, it is clear that the additional freedom is 
beneficial, since the lift constraint is now small and determinate. The constraints on 
pitching and rolling moment, which are further reduced, largely arise from small drag 
and side force constraints, which produce appreciable moment constraints at the balance 
virtual centre, about which the moment changes are measured. As pressurised zeros are 
taken before and after each run, there is no appreciable overall reduction in balance 
accuracy. Scheme B requires the block to be stable to small displacements. Except 
as far as lateral translational displacements are concerned, the system is inherently 
stable, since the local pad pressure increases automatically as the gap decreases 
at a given internal pressure, so that restoring forces and moments are produced. The 
lateral movement of the block, however, had to be restricted by a set of three balls 
running freely in sockets. 


Scheme B has now functioned successfully for more than two years. 
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Appendix III 
A Jet-PaTH VISUALISATION TECHNIQUE 


The path of deflected plane jets or circular jets can be visualised quite simply by 
means of fluid injection method. Streams of a suitable liquid, for example, water or 
paraffin, fed by gravity from a static tank reservoir above the tunnel, are injected into 
the centre of the jet just outside the slot. The liquid is dispersed as a cloud 
of fine drops, and rapidly spreads across the jet as it is carried downstream. Elevation 
or plan views of the path of the jet can be obtained with conventional lighting, and 
transverse sections by projecting a plane of light normal to the main stream direction. 


The method is particularly suitable for jet speeds above 500 ft./sec. The effect 
of model configuration changes can be studied in as much detail as required, and no 
particular precautions are necessary, apart from mopping up the liquid puddles which 
form on the tunnel floor! To visualise a flow of 1 Ib./sec. of air, from one wing 
of a jet-flap complete model with a slot area of about 300-05 in., eleven tubes of 
0-05 in. internal diameter were regularly spaced along the span, together with an extra 
tube of 0-10 in. diameter to visualise the tip vortex. About one gallon (10 Ib.) of water 
sufficed for 10 minutes’ continuous running. 
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Space probe reaches 
heights of over 500 miles— 
speeds of over Mach 10— 

with unprecedented reliability... 


...AND BRISTOL SIDDELEY 
SUPPLY THE POWER 


One of the largest manufacturers of motive 
power units in the world, Bristol Siddeley 
Engines Limited produce the Gamma. A liquid 
propellent rocket engine, the Gamma powers 
the Saunders-Roe Black Knight, Britain’s 
highly successful space research vehicle. An 
extremely reliable powerplant, the Gamma pro- 
duces a total sea-level thrust of 16,400 Ib (7,438 
kg) and nearly 19,000 Ib (8,618 kg) outside 
the earth's atmosphere, for a total powerplant 
weight of under 700 Ib. 


At the Woomera rocket range in Australia, 
the Gamma has sent Black Knight over 500 
miles into space at speeds in excess of Mach 10 
with unprecedented reliability. For, to date, 
every Black Knight vehicle has been launched 
successfully. 


Since Bristol Siddeley’s rocket division 
began work in 1946, it has developed a wide 
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range of components. By combining these com- 
ponents in single or multi-chamber layouts, 
thrust requirements from 500 lb up to 100,000 lb 
can be met. 


Bristol Siddeley products currently in production 
tnelude 
AERO: 
TURBOJETS: Olympus; Sapphire; Orpheus; 
Viper. 
TURBOPROPS: Proteus; Double Mamba. 
ROCKETS: Gamma and others. 
RAMJETS: Thor | & 2. 


INDUSTRIAL & MARINE: 
DIESELS: Maybach. 
GAS TURBINES: Proteus. 
OTHER PRODUCTS: Gas bearing compressors; 
Beaver ball screws and splines. 
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BRIsToL Aircraft 


OPPORTUNITIES FOR 
WORK ON NEW 
PROJECTS AT BRISTOL 


Bristol Aircraft Limited offer excellent and progressive 
careers with their research, development and design 
teams at FILTON to young men or women who have a 
degree or a Higher National Certificate or Diploma. 


Bristol are known throughout the world for their suc- 
cessful turboprop BRITANNIA and the BLOOD- 
HOUND Guided Missile. 

They can offer interesting and challenging careers on 
the projects of the future. 


Structures Design Engineers 
Mechanical Engineers 
Dynamics Engineers 
Mathematicians Draughtsmen 
Physicists Designers 
Aerodynamicists 
Systems Engineers 
Electrical/Electronic Engineers 


If you are interested, write to:- 

J. Raimes, Personnel Manager, 

Room FW/10/AQ 

Bristol Aircraft Limited, FILTON, Bristol, 


giving full details of your education, experience and position 
in which you are interested. 


MEMBER COMPANY OF BRITISH AIRCRAFT CORPORATION 
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if you are an A.F.R.Ae.S. or possess an equivalent 
qualification, please write to the Technical Staff Manager, 
when we shall be pleased to give you details of 
the opportunities available to experienced men within our 
technical organisation. 


Blackburn Aircraft Limited, Brough, Yorkshire 


Alzs/d 
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